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“The Physics 


TIMELY and interesting report of an after- 
noon lecture delivered by Professor J. D. 
Bernal appeared recently in ‘‘Nature.”* This 
lecture delivered at the Royal Institution was 
entitled ‘“The Physics of Air Raids’’ and has as 
its purpose the quieting of fears by a rational 
explanation of what can and what does occur 
in bomb explosions. Unknowns are always much 
more fearful than knowns. For instance, a little 
understanding of electric discharges goes a long 
way in quieting one’s nerves in a thunderstorm; 
particularly if one happens to be in a closed 
automobile. Another very important point is 
that rational means of minimizing the effects of 
air raids require first of all an understanding of 
the physics involved. 

Professor Bernal has the faculty of eliminating 
from his discussion many of the unnecessary 
details and thus making the general principles 
stand out clearly. He opens his article with an 
explanation of the fact that as a bomb explodes 
a single wave pulse or sound of very high in- 
tensity is generated. If sounds of all frequencies 
and intensities traveled at the same speed, this 
high intensity pulse would retain its shape as it is 
propagated outwards. However, since the elas- 
ticity of the medium increases considerably 
under the enormous compression involved in 
explosions the high pressure phase breaks away 
from the center of the pulse just as the tops of 
water waves break ahead near the shore. Thus 
the sharp wave of the explosion is smoothed out 
at relatively short distances from the source. 

Many seemingly freakish effects occur during 


* J. D. Bernal, ‘‘The Physics of Air Raids,’ Nature 147, 
594 (May 17, 1941). 


of Air Raids” 


air raids which are in reality fine examples of the 
application of the most elementary laws of 
physics. In a narrow street with tall buildings 
on either side, an air compression rises to. the 
top of the buildings where free expansion can 
take place. A suction wave then travels back 
toward the street. The suction wave may pull 
windows and doors out of place in a violent and 
seemingly impossible way. Another freakish 
effect is the sharpness of the shadows produced 
by the high frequency components of the intense 
wave pulse. For example, a person behind a 
garden wall close to an explosion is uninjured 
while some one else at a considerable distance 
is badly hurt. 

Since much of the damage from air raids is 
produced by splinters and earth shock, Professor 
Bernal has analyzed in a broad way the mecha- 
nism of the explosions. When the bombs burst, 
the expanding gases are held in a more or less 
spherical compression chamber. Failure takes 
place along a conical surface and fragments 
varying from 40 lb. lumps to dust particles are 
projected with speeds up to 4000 ft. per sec. 
By studying the penetration of thesé particles, 
more efficient ways of stopping them can be 
developed and the bomb shelters improved. 

There is need for more discussions of the type 
that Professor Bernal has written. Naturally the 
requirement of secrecy in connection with war 
work limits the range of topics which can be 
discussed. However, a statement of general 
principles without giving the details of the 
research valuable to the enemy serves the very 
useful purpose of quieting the nerves of the 
general public and giving the people more 
confidence in their leaders. 
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Influence of Certain Variables on the 
Stresses in Gear Teeth* 


By THOMAS J. DOLAN 


University of Illinois, Urbana, Illinois 


A photoelastic study was made to obtain a more accurate picture of the localized 
stresses in gear tooth fillets than can be obtained from the commonly used Lewis 
equation. A set of approximate stress concentration factors are presented for fillets 
in gear teeth based upon a set of variables such as the radius r of the fillets, thick- 
ness t of tooth, height h of load position on the tooth, etc. The problem was attacked 
by making a series of tests of ‘conventionalized gear teeth’’ in the form of short 
cantilever beams in which one factor at a time could be varied. The values of stress 
concentration factor were computed as the ratio of the maximum stress in the fillet to 
the calculated stress at the ‘‘assumed weakest section,” as defined by Lewis’ equa- 
tion. In general, it was noted that as the height of load position (h/t) was increased, 
the maximum stress increased, but the stress concentration factor decreased. An 
approximate equation for the tensile stress concentration factors, k, for the con- 
ventionalized models with sides of the tooth parallel is: k=1.25(t'/r)®*(t/h)®*. 
A definite decrease in the tensile stress concentration factor was observed as the 
tooth model was gradually tapered inward to approach Lewis’ ideal shape of the 
parabola of uniform strength, whereas no change in the compressive stress concen- 
tration factor occurred as the angle was varied over a 25° range. Tests now in 
progress on generated gear teeth having a 14}-degree pressure angle and a diametral 
pitch of 2 indicate that the stress concentration factor for the tensile fillet is given 


fairly accurately by the formula: K =0.22(t/r)%*(t/h)**. 


Introduction 


OCALIZED stresses of large magnitude are 
always present in a loaded gear tooth at the 
point of load application and at the fillets. In 
order to facilitate and improve present methods 
of gear design there is a definite need for more 
exact information as to the magnitudes of these 
stress concentrations. The present paper will be 
limited to a discussion of some of the factors 
that affect the stress concentrations occurring at 
the fillets of a gear tooth, as determined by a 
series of photoelastic tests on models designed to 
‘simulate the structural action of spur gear teeth. 
Investigations have been carried out along some- 
what similar lines by Timoshenko,' Baud and 
* Paper presented at the Eleventh Semi-Annual Eastern 
Photoelasticity Conference, May 24-25, 1940, at the Car- 
negie Institute of Technology, Pittsburgh, Pennsylvania. 


1S. Timoshenko and R. V. Baud, “Strength of gear 
teeth,”” Mech. Eng. 48, 1105 (1926). 
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Peterson,*? Black,‘ and others; it was hoped 
that the present investigation might cover a 
wider range of variables and yield a clearer 
understanding of the relative effect of some of 
the variables affecting the localized stresses at 
gear tooth fillets than the results thus far avail- 
able. It is realized, of course, that the photo- 
elastic method yields only theoretical values of 
stress based on elastic action in the members; no 
attempt has been made to determine the extent 
to which these theoretical stresses are of sig- 
nificance in causing damage to actual metal 
gears in which the relative ‘‘notch sensitivity,” 
localized plastic readjustments in the material, 
and distribution of the loads among several 





2 R. V. Baud and E. Hall, ‘‘Stress cycles in gear teeth,” 
Mech. Eng. 53, 207 (1931). 

3R. V. Baud and R. E. Peterson, ‘“‘Load and stress 
cycles in gear teeth,’’ Mech. Eng. 51, 653 (1929). 

*P. H. Black, ‘‘An investigation of relative stresses in 
solid spur gears by the photoelastic method,” Bull. 228, 
Eng. Exp. Sta., University of Illinois (1936). 


JOURNAL OF APPLIED PHYSICS 














teeth may greatly influence the damaging value 
of stress. 

At present most -designers use the Lewis 
equation, in designing spur gear teeth, or employ 
modifications of this equation in the design of 
bevel gear, helical gear, and other special forms 
of gear teeth. The basic assumptions as made by 
Lewis in the equation for the maximum stress 
developed at the fillet in a gear tooth were: 
(a) the load is carried by one tooth and acts 
along the face of the gear at the tip of the tooth 
in a direction normal to the tooth surface; 
(b) this force can be resolved into two com- 
ponents at the intersection of the line of force 
with the radial center line of the tooth, one com- 
ponent acting radially and the other tangentially ; 
the radial component was disregarded in his 
analysis because it was small and served as (an 
unreliable) source of strength to materials weak 
in tension. The Lewis equation treats the gear 
tooth as a cantilever beam and determines the 
value of a nominal flexural stress at an ‘“‘assumed 
weakest section’’ of the tooth. This “‘weakest 
section”’ is located between the two points of 
tangency of the tooth profile with a parabola 
inscribing in the tooth outline (referred to as the 
parabola of uniform strength) with its apex at 
the intersection of the action line of the applied 
force and the radial center line of the tooth. 
For example, in Fig. 1, the parabola ACB locates 
the “weakest section’”’ between the points of 
tangency A and B. Since it is based on erroneous 
assumptions, this equation does not give the true 
stress values developed, but has served remark- 
ably well for a period of years without change 
(probably due largely to the accommodating 
properties of the materials used and to the low 
working stresses employed in actual design). 


Object and Scope of the Tests 


A photoelastic study was made in an attempt 
to obtain a more accurate picture of the localized 
stresses in gear tooth fillets than can be obtained 
from the Lewis equation. The main purpose of 
the tests was to determine a set of approximate 
stress concentration factors for fillets in gear 
teeth based upon a set of variables such as the 


5 QO. A. Leutwiler, Elements of Machine Design (McGraw- 
Hill Book Company, Inc., New York, 1917), Chapter XII. 
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radius r of the fillets, thickness ¢ of tooth, height 
h of load position on the tooth, and the angles 
at which the flank and bottom land of the tooth 
meet; these variables were selected as being the 
ones that were significant in determining the 
effect of the shape of the tooth on the localized 
stresses developed. 


C= 





















































Fic. 1. Gear tooth notation and model dimensions 
A and B are points of tangency of parabola ACB with the 
tooth profile. 


In order to present the results in a form most 
familiar to the gear designer, it was found ad- 
visable to express the localized stresses in the 
gear teeth (or the equivalent test specimens) in 
terms of the Lewis equation but modified by the 
introduction of a stress concentration factor. 

Before proceeding with an analysis of the 
stresses in generated gear teeth, it was deemed 
advisable to study the effect of the shape factors 
that are altered in producing the various con- 
tours of generated gear teeth to determine their 
relative importance in influencing the magnitude 
of the stress concentration factor. This phase of 
the problem was attacked by making a series of 
tests of ‘“‘conventionalized gear teeth’’ in the 
form of short cantilever beams in which one 
factor at a time could be varied. 

The series of alterations of shape that were 
chosen to represent the most important variables 
for this series of tests of ‘“‘equivalent’’ models 
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are illustrated diagrammatically in Fig. 1. These 
may be briefly summarized as follows: 


(1) The radius rof the fillet between the flank of the tooth 
and the root circle. 

(2) The nominal thickness ¢’, and the thickness ¢, at the 
‘‘assumed weakest section”’ of the tooth. 

(3) The height of load position 4 above the “assumed 
weakest section.” 

(4) The angle 8 between the side of the tooth and the 
radial centerline. 

(5) The angle a between the bottom land of tooth space 
and the thickness chord. 

(6) The angle ¢ that the action line of the load makes 
with the tooth face. (This angle was maintained 
constant at 90° in the tests.) 





A series of models was made with the dimen- 
sions shown in Fig. 1 having nominal values of 
L=3", D=3", and U =1”", and values of h were 
varied from }” to 13”. All models were made 
from {” thick polished plates of BT61893 Bakelite 
and were annealed after roughly sawing them 
slightly oversize from the plate. Immediately 
after annealing, they were machined to finished 
size with light cuts, and the tests completed 
within a few hours after machining. No appreci- 
able residual stresses or edge effects were present 
in the models tested. Calibration beams used for 
determination of the fringe stress values were 
cut from the rear portion 
of the model after the tests 

















Fic. 2. Typical fringe photographs of several models. 
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on the model were com- 
pleted. 

Each model was _ sup- 
ported by lightly clamping 
(the lower edge of length 
“L” in Fig. 1) to a rigid 
bar, and was loaded with a 
line load across the face of 
the tooth ; successive fringe 
photographs were made in 
monochromatic (mercury 
green line) light with a 
series of different load 
positions h. The model 
was then altered by vary- 
ing the fillet radius r (or 
the angle 8 of the tooth 
sides, etc.,) in a series of 
successive increments, and 





a complete new set of fringe 
photographs were made 
with loads at the same 
positions for each new al- 
teration of the model. Eight 
of these ‘‘equivalent’’ mod- 
els were tested, and approx- 
imately 170 negatives of 
fringe photographs were 
obtained and analyzed. 


Test Data and Discus- 
sion of Results 








In Fig. 2 are shown 
photographs of typical 
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MODEL 6 
r= 0.25) t=1.09", r/t'=0.229, h/t=0.557 


Fic. 3. Stress distribution curves. Stress scale: 1’/°=10 
p.s.i. Approximate enlargement =0.9/1. Curves based on 
l-inch face. 


fringe patterns for several of the models tested, 
and in Fig. 3 and Fig. 4 are shown typical dis- 
tributions of stress around the fillets for several 
of the test loadings. In analyzing the results, 
the stresses were converted to values of stress 
for a ‘unit’? model (1” thick) loaded with a 
unit load (one pound), and the nominal values 
for this same ‘‘unit’’. model were computed from 
the Lewis equation. It is a relatively simple 
matter to obtain the value of this calculated 
stress for any load position by the use of a 
graphical method of locating the assumed weak- 
est section, and a direct application of the 
flexure formula. One set of results showing 
changes produced in the fillet stresses by a 
variation of the load position and fillet radius 
is plotted in Fig. 5, in which all other variables 
of shape were maintained constant. 

It will be observed that the stresses increased 
rapidly as the fillet radius was reduced, and that 
the stresses also increased (but not in direct 
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proportion) as the height of load position above 
the bottom of the tooth was increased. Further- 
more, the maximum stress on the fillet in com- 
pression was, in general, often larger than that 
at the fillet in tension. From the viewpoint of 
design, the maximum tensile stress is perhaps the 
most important consideration, but the stress 
concentration factors for both fillets of each 
tooth were computed and these values are shown 
in Fig. 6 for the above data. 

The values of the stress concentration factors 
were computed as the ratio of the actual maximum 
stress in the fillet, as obtained from the photoelastic 
fringe pattern, to the calculated stress at the 
“assumed weakest section,”’ as defined by Lewis, 
using his original assumption that the radial com- 
ponent of load can be neglected. 

In general, it will be noted that as the height 
of load position (h/t) was increased, the maxi- 
mum stress increased, but the stress concentration 
factor decreased. It is also evident that a generous 
radius on the fillet was very effective in de- 
creasing the stress concentration factor. 

By scaling off the ordinates to the curves in 





MODEL 6 
r= 0.25) t=1.09) r/'=0.229, h/t=0557 








~5.68 


MODEL 7 
r= 0.28" 1'=1.258" r/t'= 0.223, h/t: 0545 
Fic. 4. Stress distribution curves. Stress scale: 1’’=10 


p.s.i. Curves based on 1-inch face. Approximate enlarge- 
ment =0.78/1. 
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Fic. 5. Variation of the fillet stresses with the tooth fillet 
radius and the load position. 


Fig. 6 at various ratios of h/t, the data may be 
replotted as shown in Fig. 7; for comparison 
with these results, curves determined by Frocht® 
and by Weibel’ have been added showing stress 
concentration factors for a beam with fillets and 
subjected to pure bending. This type of loading 
would approach the condition for h/t= ~ and it 
will be noted that the localized stresses were 
much more severe in the short stubby cantilever 
beams that approximated the proportions of 
gear teeth. An approximate equation for the 
tensile stress concentration factors k for the 
conventionalized models with sides of the tooth 
parallel is: 
k=1.25(t'/r)%*(t/h)*?. 

*N. C. Riggs and M. M. Frocht, Strength of Materials 

(Ronold Press Company, New York, 1938), p. 389. 


7E. E. Weibel, ‘Studies in photoelastic stress determina- 
tion,” Trans. A. S. M. E. 56, 637 (1934). 
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This relationship fits the data within an accuracy 
of about 5 percent for the range of values covered 
in the tests. 

The curves in Figs. 8 to 10 show the variations 
obtained for the stress concentration factor when 
the fillet radius was kept constant but the angle 8 
of the side of the tooth varied. A definite de- 
crease in the /ensile stress concentration factor 
was observed as the tooth model was gradually 
tapered inward to approach Lewis’ ideal shape 
of the parabola of uniform strength, whereas no 
change in the compressive stress concentration 
factor occurred as the angle was varied over a 
25° range. 

Tests of models in which the angle a, between 
the bottom land of tooth space and the thickness 
chord, was varied, showed no measurable change 
in stress concentration at either the tensile or 
the compressive fillet for angles @ up to 15 
degrees. This was perhaps to be expected since 
the material removed from the model by in- 
creasing the angle @ was usually in regions of 
relatively low stress and consequently did not 
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Fic. 6. Variation of the stress concentration factors with 
the tooth fillet radius and the load position. 
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“1G. 7(a). Variation of the tensile stress concentration 
factor with the load position and the tooth fillet radius. 
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Fic. 8. Variation of the tensile stress concentration factor 
with the angle 8 and the load position. 
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Fic. 7(b). Variation of the compressive stress concentra- 
tion factor with the load position and the tooth fillet 
radius. 
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Fic. 9. Variation of the tensile stress concentration factor 
with the load position and the angle 8. 


necessitate much readjustment of stresses in the 
regions of high stress. 

A companion series of tests is now in progress 
in which models of actual generated gear tooth 
shapes are being studied. However, the variations 
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in shape of the models in this latter series of 
tests are being controlled by the process of 
generation of the tooth to conform to standard 
gear practice. For instance, in Fig. 11 the tooth 
is generated by a rack tool rolling on the pitch 
circle, and hence the shape depends upon the 
pitch, edge radius of the cutting tool, addendum 
and dedendum lengths, and the pressure angle. 

In comparing the preliminary results of tests 
on actual tooth forms with the data shown in 
Fig. 7, one important variation between the two 
sets of tests has been noted, namely: the tensile 
stress concentration factors (shown in Fig. 7(a)), 
were from 12 to 25 percent higher than those 
obtained for models of true gear tooth contours 
having the same minimum fillet radius and 
relative proportions, whereas the compressive 
stress concentration factors (shown in Fig. 7(b)) 
ranged from zero to ten percent higher than those 
obtained from models of actual gear teeth. The 
greatest differences between the two sets of data 
occurred with the fillets of sharpest curvature 
(small r/t/ values). 

One reason the stress concentration factor for 
the tensile fillet in the actual gear tooth was 
much smaller than the values shown in Fig. 7(a) 
is that the generated gear tooth models were 
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Fic. 10. Variation of the compressive stress concentration 
factor with the angle 8 and the load position. 
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GEAR TOOTH NOTATION 
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fa | 

\ -t j—+|' / 
| GEAR TOOTH } 
en. ee 


Fic. 11. Gear tooth notation. Angle g= 14}°. 





loaded in a direction normal to the tooth face, 
and the radial (compressive) component of load 
was neglected in the analysis by using as the 
theoretical stress only the value of Lewis’ flexural 
stress. This small radial component of load would 
produce a slight lowering of the observed stress 
at the fillet in tension, and hence result in a 
smaller stress concentration factor. Similarly, the 
radial component of load would produce a higher 
observed stress at the fillet in compression, and 
should result in a higher stress concentration 
factor being obtained. However, a rough check 
on the radial components of load indicated that 
their magnitudes were small and they produced 
only a very small increment of stress as com- 
pared with the flexural stress developed by the 
transverse component of the applied load. 

The main reason that the generated gear teeth 
exhibited somewhat lower stress concentration 
factors than did the conventionalized models 
probably lies in the fact that the conventional- 
ized tooth models had circular fillets, whereas the 
fillets on the true gear tooth were somewhat 
spiralled by the process of generation with the 
minimum radius at the root being the assumed 
value of the fillet radius. Consequently, it is 
evident that the equivalent circular radius of the 
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fillet on a generated gear tooth that is effective 
as a stress raiser is not the minimum radius of 
curvature of the fillet contour; if the values of 
stress concentration factor shown in this paper 
were used for design purposes they would prob- 
ably be on the safe side and yield values of stress 
that were slightly too high. 

Further tests are now in progress on the 
stresses in models of actual generated gear teeth, 
covering a wide range of standard shapes. It is 
hoped that, by a careful analysis of the results, 
the above variations can be accounted for and a 
set of stress concentration factors, based on the 
shape factors of the teeth, can be devised which 
will be an aid to more efficient gear design. The 
tentative results of the tests thus far on generated 
gear teeth having a 143-degree pressure angle, 
and a diametral pitch of 2, indicate that the 
stress concentration factor for the tensile fillet as 
obtained by dividing the observed stress by the 
computed combined (flexural plus direct) stress 
is given fairly accurately by the formula: 


K =0.22(t/r)%?(t/h)* 4. 


In closing, it should again be emphasized that 
the significance of photoelastically determined 
localized stresses for purposes of design of actual 
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metal members is greatly influenced by the 
properties of the materials of which the members 
are made; the exact importance of these theo- 
retical stresses in measuring the nearness to 
structural damage of the member can only be 
determined by correlation with carefully con- 
ducted laboratory fatigue tests and by com- 
parison with service records. 
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On Similarities Between Stress and Flow Patterns * 
By MIKLOS HETENYI 
Westinghouse Research Laboratories, East Pittsburgh, Pennsylvania 
(Received March 29, 1941) 


The purpose of this paper is to review established relationships between stress 
and flow problems. Two-dimensional plane stress systems were found to have 
features analogous to non-viscous potential flow fields on one hand, and to viscous 
flow on the other hand. Accordingly, this discussion is divided into two parts. 
In each case the fluid is assumed to be continuous, homogeneous, and in- 


compressible. 


I. Similarities Between Plane Stress and 


Potential Flow Fields 


i. In a plane stress system there exist the fol- 
lowing relations between the stress components 
Tx, Oy, Try, the strain components e,, €,, Y., and 
the displacement components u and v: 


Ou 1 

é,= =—(o,—wG7,), 
Ox &E 
Ovi 

€, = — (o,—po,), 
oy ke 
Ou. dv 2(1+yp) 

Y: uu _ + = T; ue 


Oy Ox Dp 


where E denotes the modulus of elasticity and p 
the Poisson’s ratio of the material. Substituting 
these values into the equilibrium equations: 


00,/OX+0T,,/ 0V=90, 


00, OV+0r,, OX=0, 
we have: 


2 dfodu_ dv 0 sd0v Ou 
1—yp 0x Ox oy Oy\0x dy 
2 d/f0u dv 0 fdv Ou 
+~)+—(—-—) -0 
1—pdy\dx dy Ox\Ox dy | 
Putting 
Ou/dx+0v0/dy=(1—p)A/ pm, 
0v/ dx —du/dy= 2a, 


* Paper presented at the Eleventh Semi-Annual Eastern 
Photoelasticity Conference, May 24-25, 1940, at the Car- 


negie Institute of Technology, Pittsburgh, Pennsylvania. 
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where \ denotes the lateral contraction and w 
the rotation at any point of the stress field, it is 
seen that 


O*vX Or 0-w JO-w 
+——=0; +—=0 
Ox* dy Ox* dy" 
while . a re 
OX 0 OX 0 
=—w; = ——w 
Ox p Oy OV bu Ox 


Comparing these results with a two-dimensional 
flow problem where the velocities v,= 0¢/dx and 
v, =0/dy can be considered as derivatives of a 
velocity potential ¢, the conditions of a source- 
free and irrotational 
dv,/dx+dv,/dy=0 and 
spectively. Since the relations between the stream 
function y and the velocity potential ¢ are 


written as 
0v,/ 0x — 0v,/dyV=0, re- 


flow can be 


Om Ox=d0Y~/0vV; I6/dvV= —dOY, Ox 
in addition to €°6/dx°+0°o/dy? =0 ,. (2) 
we also have 0°y/dx°+07*~/dy" =0 


Comparing Eq. (2) with Eq. (1) it is seen that 
the relationship between ¢ and y is the same as 
that found between \/yu and w. Hence the con- 
clusion that: 

The equipotential and stream lines of a flow 
pattern can always be considered as lines of equal 
lateral contractions (isopachics) and lines of equal 
rotations in a corresponding plane stress system. 

Possibilities of using this analogy in solving 
stress problems were discussed by R. Doerfel.' 

2. Another analogy has been established by 
relating elastic displacements to flow velocities. 
Using the above-introduced notations and taking 
the displacement components u and v propor- 
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tional to the velocity components v, and v, 


u=kv,=kdy/dy and v=kv,= —kdy/dx, 


we find by substitution that these expressions 
will satisfy Eq. (A). Hence: 

The velocities of a potential flow can be con- 
sidered as displacements in a plane stress system. 

It is to be noted that the opposite is not true. 
From the fact that in every stress system equiva- 
lent to a flow pattern the displacements will 
have to be simple harmonic functions, it follows 
that such stress systems where the displacement 
components are not simple harmonic functions 
will have no equivalents in potential flow. 

This displacement-velocity analogy is dis- 
cussed in a paper by J. Barta? with reference to 
aeronautical problems. 

3. The conditions under which the orthogonal 
system of the equipotential and stream lines in a 
flow problem is equivalent to a system of stress 
trajectories were investigated by P. Neményi® 
and U. Wegner.‘ Their findings can be sum- 
marized as follows: 

A. It is always possible to establish a system of 
stress trajectories equivalent to the stream and 
potential lines of a given flow problem. 

Denoting by w(z) the known complex velocity 
function (s=x+iy) which is the derivative of a 
complex potential w(z)=S’(s) the equivalent 
stress system can be found as 


o,=7Re[w*(s) |— 2rd f (at-eaar(at+p], 


zu 


o,=— nRe{o%(2)]-2R¢ f (atta)ar(al+p], 


Try = — nl m{[w*(s) |. 


Here Re and Jm denote the real and the imagi- 
nary part of the complex function in question 
and »=a(x°+ y")+2bx+2cv+d; a=b+ic. There 
are five real parameters in the above formulas: 
a, b, c, d and p. Consequently there is a variety 
of stress systems whose trajectories are equiva- 
lent to those of the given flow problem. The 
physical interpretation of the parameters is as 
follows: a, b and c are geometric constants; a de- 
termines the scale, b and c determine the origin 
of the coordinates. The parameter d shows the 
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possibility of superposing a harmonic stress field 
which does not change the system of stress 
trajectories while p gives a superposed uniform 
tension or compression. 

B. The solution for finding a flow problem whose 
stream and potential lines are equivalent to the 
trajectories of a given stress system is possible only 
if one of the following conditions is satisfied: 

2(az+d) 
(1) N(s)= F(z), 
az+a 


0” 0" 
(2) ( + ) toe o1—o2 
Ox* dy? 


4(ad — b? —c?) 
~ Ca(x?-+-y2) +2bx-+2cy-+d}? 
F(z) | [a(x*?+-y?) 
+2bx+2cy+d], 
o,=+Re[ N(z)], 
try =Im[N(z) }. 


ey" esi > 


o,= —Re[ N(z)], 


Here &=b—ic and N(z) and F(z) are two ana- 
lytical functions of the complex variable , 
whose imaginary parts have to satisfy the equa- 
tion 


Try = (x?+y?) Im[ F(z) ]+ImLN(z) J. 


Otherwise the notations are the same as above. 
Wegner! suggested using condition (2) in con- 
nections with photoelastic test results. The pro- 
cedure would then be to determine photoelasti- 
cally at every point of the stress field the value 

o,—o2, then carry out the assigned operations 
on the left side of Eq. (2), and finally to try out 
if, by a proper choice of the parameters a, b, ¢ 
and d on the right side, the equality could be 
satistied. This semi-experimental method seems 
that its final result could 
hardly be accepted for an accurate proof. 

In both of the above cases, A and B, the ques- 
tion of boundary conditions was still left open. 
In general, we will find that the boundaries of 
the equivalent stress field will have to be sub- 
jected to distributed loadings of varying in- 
tensity. Of greater interest are the cases where 
an equivalent stress system with unloaded 
boundaries can be found. This phase of the 


to be so involved 
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problem was investigated by H. Neuber® who 
also disclosed some new and interesting properties 
of the stress trajectories. He has shown that, if 
the orthogonal trajectories of the principal 
stresses ¢, and o2 are given as two sets of curves, 
p=const. and g=const., where p=/p(x, y) and 
qg=q(x, y) are functions of the Cartesian coordi- 
nates x and y, then the equilibrium equations can 
be written as 


h,da,/dqg+(o1—02)dh,/dq=0\ 


? (i 
h,002/0p —(a,—02)dh,/dp= Of v 


Here h, and h, denote the spacing of the curves 
p and q, respectively. Assuming that a curve 
p=const. represents an unloaded boundary we 
put o.=0 and find that, in this case, the first 
of the above equations reduces to 


h,d0,/0q+0,0h,/dqg=9/0q(o,h,) =0, 
the solution of which is 
o;1> . h p* 


According to this result, which is of general 
validity and is not limited to analogous stress 
and flow fields: In any plane stress system the 
variable distance h, between the unloaded boundary 
and the principal stress trajectory p is inversely 
proportional to the boundary stress value o;. 

It is seen that on the basis of this fundamental 
relationship alone a new optical method of 
stress analysis could be developed, which would 
make use only of the isoclinic lines, deriving 
from those the stress trajectories and obtaining 
thus the principal stress values and stress con- 
centration factors at any point along the un- 
loaded boundary. 

Many of the problems related to stress and 
flow fields can be analyzed most conveniently in 
the curvilinear coordinate system p and gq. We 
find, for instance, that if the principal stress 
trajectories p and q are to be considered as stream 
and potential functions of an analogous flow 
‘problem, then they have to satisfy the Cauchy- 
Riemann equations 


0q/dx=dp/dy; dq/dy=—dAp/dx 
and also 
dx/dqg=dy/dp; dx/dp=—dy/dq. 
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A substitution of these expressions into the 
formulas 


h,? = (dx /dq)*+ (dy /dq)* 
and h,*> =(dx/dp)*+(dy/dp)* 
leads to 

h,»=h =h. 


An orthogonal system of trajectories which has 
the above property is called isometric. Hence we 
have the general rule that: 

Only an isometric set of stress trajectories can be 


considered as a set of potential and stream lines of 


a flow problem. 
From the conditions of equilibrium Eq. (a) 
and the compatibility equation 


A(oi+02) =() 
follows by integration that 
o,—02=(1/h?)[a(x*+y*) +bx+cyt+d ], 


which, along the unloaded boundary where o2=0, 
through the substitution ¢;=C/h gives the for- 
mula 


o,=C?/(a(x?+y") +bx+cy+d ]. 


Thus, in an isometric case, when the stress trajec- 
tories can be considered as potential and stream 
lines of a flow pattern, the principal stress along the 
unloaded boundary must be expressible in the above 
form. This formula can be used for deciding 
whether a given stress problem with known 
boundary stress values has an equivalent flow 
field. 

From the above formula we also find that the 
principal stress along the unloaded boundary will 
vanish if C=0, and then it can be shown that also 


a(x*+y") +bx+cy+d=0. 


Since the above equation is the general equation 
of a circle, the conclusion is that: 

In an isometric case the stress along the unloaded 
boundary will vanish only if the boundary is a 
circular segment or, in the limiting case, a straight 
line. 

Neuber has shown a variety of stress problems 
which are isometric and thus have an equivalent 
flow pattern. The boundary lines in these 
problems are composed of circles, which are 
first assumed to intersect giving various crescent- 
shaped figures, and secondly are assumed not to 
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intersect and result in plates having circular 
holes. As a special case of the first class of these 
problems we have the diametrically compressed 
circular roller, where both of the principal 
stresses vanish along the boundary and which 
has been recognized for long as a case which 
has an equivalent flow pattern. 

4. A number of analogies were established by 
Lord Kelvin, Boussinesq and others® between 
flow problems of various types and the stress 
distribution in a prismatical bar subjected to 
torsion. Though problems of torsion do not 
belong to the class of plane stress, they may be 
mentioned here on account of their two-dimen- 
sional character. On the basis of the Kelvin 
analogy, experiments were carried out by Den 
Hartog and McGivern.? 


II. Similarities Between Plane Stress and 
Viscous Flow Problems 


The stream function of flow ¥, whose deriva- 
tives give the velocity components v,= —dy/dy 
and v,=dy/dx, satisfies the equation AAy=0 
where A=0?/0x*+0?/dy? is the Laplacian oper- 
ator. The flow of the viscous fluid is assumed 
infinitely slow. The same biharmonic equation is 
satisfied by the transverse deflection of an 
elastic plate, and also by the Airy stress function 
of any plane stress system. The first of these 
analogies was established by Lord Rayleigh® 
while the second one, which will be discussed 
here, has been applied by J. N. Goodier.® 


o:=0°F/dy*, o,=0°F/dx*, 12y= —0°F/dxdy, 


we have 
AF=o,+0,, 


which is proportional to the lateral contraction 
in a plane stress problem. In case of viscous flow 


Ay = dv, /0x —dv,/dy 


represents the vorticity of the flow. Hence we 
have the general correlations that: lines of 
constant thicknesses (isopachics) in the stress 
problem will correspond to lines of constant 
vorticity in the flow field and the lines of con- 
stant rotations, being orthogonal trajectories of 
the isopachics, will be equivalent to the lines of 
constant fluid pressure in the flow problem. It 
also follows that the stress trajectories and iso- 
chromatics will represent in the field of flow 
the orientation of maximum distortion and the 
lines of constant maximum distortion, respec- 
tively. 

In this analogy the equivalent boundary con- 
ditions have also simple physical interpretation. 
If the rigid boundaries in the flow problem are at 
rest or in uniform motion the corresponding 
stress system will have the same boundaries 
stress free or subjected only to uniform normal 
stress. The tangential drag and normal pressure 
along the boundaries of the flow field can be 
determined from the isochromatic and isopachic 
lines of a corresponding photoelastic model. 
Several interesting applications of this analogy 


Since the second derivatives of the stress have been suggested by Goodier in his above- 
function F give the stress components mentioned paper. 
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A Rational Approach to the Numerical Solution 
of Laplace’s Equation * 
( The Four-Point Influence Method ) 


BY M. M. FROCHT** AND M. M. LEVENT 


Carnegie Institute of Technology, Pittsburgh, Pennsylvania 


HE solution of many physical problems 

depends upon the numerical solution of 
Laplace’s equation for given boundary values. 
Methods have been perfected which materially 
reduce the time and labor involved in the itera- 
tive process. In this process arbitrary values are 
assumed at the points of intersection of a square 
network, and these values are systematically and 
repeatedly adjusted in accordance with the 
well-known requirement ‘that the value of a 
harmonic function at every point must equal the 
mean of the values at four equidistant nearby 
points. Notable improvements in the process of 
iteration were made by Shortley and Weller.' 
They introduced the use of the nine-block unit, 
the preliminary coarse network and the general 
rectangular network. 

In this paper a method is suggested which 
further improves the iterative process by provid- 
ing a rational set of initial values in the place of 
the arbitrary numbers with which the iteration 
process generally begins. Another object of this 
paper is to provide a method for a first approxi- 
mation to the value of a harmonic function at an 
isolated interior point from given boundary 
values. Such a method is of value in problems of 
elasticity and may be likened in its purpose to a 
lateral extensometer. A simplified method is also 
given for the transitions in a mixed network. Two 
examples illustrating the method are included. 


* Paper presented at the Eleventh Semi-Annual Eastern 


Photoelasticity Conference, May 24-25, 1940, at* the 
Carnegie Institute of Technology, Pittsburgh, Pennsyl- 
vania. 

** Associate Professor of Mechanics, Carnegie Institute 
of Technology. 

+ Mechanician and part-time instructor, 
Department, Carnegie Institute of Technology. 

1See G. H. Shortley and R. Weller, ‘‘The numerical 
solution of Laplace’s equation,” J. App. Phys. 9, 334 
(1938). 
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The Four-Point Influence Equation 


Let Ya, Yo, Ye, and Yu represent the known 
values of a harmonic function (x, y) at four 
points A, B, C, and D, of coordinates (a, 0), 
(0,5), (—c, 0), and (0, —d), respectively, as 
shown in Fig. 1. It can be shown that the values 
of the function (x, y) at four such given points 
determine a first approximation to the value of 
the function at the point of intersection of the 
two lines on which the points are located, point O, 
Fig. 1, and that this approximation is given by? 


1 1 1 1 
( + ve —Wat vi 
ac bd a(a+c) b(b+d) 
1 1 


+ Vet Wu. (1) 
c(c+a) 


d(d+b) 


that the denominators on the 


right side of the equation represent simple areas. 
Thus a(a+c), c(¢+a) represent, respectively, 


a 


@ (0, b) 


We observe 


Co— 
(-Cy 0) 





=i eA —> 
0 (a, 0) 





® (0,-d) 
D 9 


Fic. 1. Representation of the known values of a harmonic 
function at four points. 


2 See reference 1, Eq. (23). 
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Fic. 2. Areas of diffusion. 
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areas A, and A, shown in Fig. 2. These areas we 
term areas of diffusion. Each of the four boundary 
points may thus be associated with an area of 
diffusion. The ratio of boundary value at a point 
by the corresponding area of diffusion we define 
as the diffused boundary value and denote by wj. 
Thus 


Wa = Wa /A ay 


etc. 
In terms of this notation Eq. (1) becomes 
Yo = B(WatwW+we+wua), (2) 


in which 
1 

(1/ac+1/bd) 
The coefficient 6 will be spoken of as the four- 
point coefficient and the value of y given by 
iq. (2) as the four-point influence. We thus have 
the following proposition : the four-point influence 
at a point is given by the product of the four-point 
coefficient and the sum of the four diffused boundary 
values. We note that the diffused boundary value 
diminishes as_ the from 


distance the origin 


increases. 


The Four-Point Influence Curve 


The four-point influence given by Eq. (1) 
represents a close approximation to the true 
value of y, at the origin O when the variation in 
the boundary about the point are 
negligible. This is, of course, always the case 


stresses 
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when the distances from the point to the bound- 
aries are infinitesimally small. When these dis- 
tances are not small the four-point influence may 
differ radically from the true local value of 
¥(x, y). If, however, the influences of several 
four-point sets from different directions be deter- 
mined, and a curve plotted with these influences 
as ordinates and angular displacements as 
abscissae as shown in Fig. 6, it is reasonable to 
assume that the mean value of these four-point 
influences will be a better approximation to the 
true value of the harmonic function (x, y) at 
the given point than any one of the separate 
values. The curve determined from the four- 
point influences will be referred to as the four- 
point influence curve, or as the influence curve. 

The influence curve repeats itself when the 
coordinate axes rotate through an angular region 
of 2/2. In the special case when a point lies on a 
line of symmetry the influence curve for that 
point becomes symmetrical at 6=2/4 where 6 
is measured from the line of symmetry. The mean 
ordinate of the influence curve furnishes the 
first approximation to the value of a harmonic 
function at an interior point from given boundary 
values. 

Example 1.—Consider the bar with the circular 
hole shown in Fig. 3, and let it be required to 


Uniform stress =3.11 fringes 
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Fic. 3. Sketch showing dimensions of a bar with an axial 
circular hole in tension and key points. 
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TABLE I. Bar with central circular hole subjected to pure tension data for key-point calculations. 


ANGULAR 
POINT POSITION a r b d a 

0° 0.392 1 3 3 9815 
14° 0.405 1 3 3 9280 
1 28.5° 0.455 1 3 1.8 7400 
37° 0.64 1 3 1.31 800 
gg 3 1 3 1.31 3270 
45° 3 1 3 1 3110 
0° 1 3 2 4 3750 
20° 1 3 2 2.75 4120 
20° 1 1.21 2 2.75 4120 
22° 1 1.00 2 2.48 4140 
2 45° 1 0.57 2 1 3500 
68.5° 2 1.06 3 1 3110 
70.5° 2 1.22 3 1 3110 
70.5° 2 4 3 1 3110 

90° 


determine the sum of the principal stresses at the 
points 1, 2, 3 and 4. It is well known that in a 
plane stress system the sum of the principal 
stresses p and q satisfies Laplace’s equation,’ i.e., 


#(p+q)/dx?+0(p+q) /dy*=0. 


We also observe that the edges of the bar are 
free boundaries, so that at every point one 
principal stress vanishes and the remaining 
stresses are parallel to the edges. Furthermore, 
the photoelastic stress pattern, Fig. 4, furnishes 
the values of (p—gq) at every point, and at free 
boundaries it gives also directly the sum of these 
stresses. 

We then proceed as follows: through point (1) 
we draw several pairs of coordinate axes, each 
pair in different directions. In the case under 
consideration these directions were chosen at 0°, 
14°, 28°, 30°, 37° and 45° with the horizontal, 
Fig. 5. The coordinate axes are produced to the 
free boundaries and the distances from the 
origin, point (1), to the points of intersection 
with the boundaries measured for each set of 
axes. These distances represent the coordinates 
a,-b, c and d entering into Eq. (1) and they, as 
well as the boundary values of (p+q), are shown 
in Table I. From these data the four-point 
influence for each set of axes is calculated by 
means of Eq. (1) and an influence curve plotted, 
Fig. 6. The area under this curve is then measured 


*For a proof see Max M. Frocht, “Isopachic stress 
patterns,” J. App. Phys. 10, 248 (1939) or Theory of 
Elasticity by S. Timoshenko (McGraw-Hill Book Com- 
pany, New York, 1934). 
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b+4q FROM 
MEAN HARMON- 
¢c b d i (p+q)i IZATION 
2890 3110 3110 7665 
2940 3110 3110 7265 
3190 3110 3750 5880 — -_ 
3400 3110 ©4100 2099 899 = 5510 
3400 3110 4100 3550 
3750 3110 3750 3590 
3750 3110 3110 3575 
2900 3110 3750 3665 
— 2120 3110 3750 1670 
— 4000 3110 4150 660 
2860 3110 2890 3064 3300 3175 
9815 3300 3280 5750 
8100 3420 3310 4920 
3110 3420 3310 3280 
3575 


with a planimeter and the mean ordinate cal- 
culated. The abrupt drops or changes in the 
curve correspond to positions of the axes at 
which there are two different boundary values. 
This happens where @=37°. One axis is then 
tangent to the circle representing the hole and 


a 





Fic. 4. Photoelastic stress pattern for the bar shown 
in Fig. 3. Model fringe value=340 psi tension or com- 
pression. 
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Fic. 5. Drawing showing coordinates and their directions F1G. 7. Drawing showing coordinates and their directions 









for key points (1) and (2). 



































for key points (3) and (4). 
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Fic. 6. Influence curves for key points (1) and (2). 
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Fic. 8. Influence curves for key points (3) and (4). 
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(p+q) equals 800 at the point of tangency or 


3270 at the right edge. At point (1) the mean 
value of (p+q) obtained from the four-point 
influence method is 5850 psi. 

In a similar manner the mean four-point 
influence value was determined for points 2, 3 
and 4 which were found to be 3300, 300 and 2740 
psi, respectively. The coordinate axes and the 
resulting curves are shown in Figs. 5, 6, 7 and 8. 

It is interesting to compare the values so ob- 
tained with what may be called theoretical 
values for the same problem which are 5510, 
3175, 465 and 2635 psi for points 1, 2, 3, and 4, 
respectively. These have been obtained by a 
process of iteration and may be looked upon as 
values at convergence. It is seen that the dif- 
ferences between the two sets of values are small. 
Moreover, the results from this problem seem to 
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Known boundary valves 
Key points 
Intervor valves obtained from normal Liebmans 


x+@0 


interwor valves obtamed from diagonal Lvedmans 


First value m each block obtamed from hey points and boundary values using 
Liebman equation with the exception of pomts C7, 08, B6, D9, and AG 
whch were calculated from the Four-pomt influence equation. 


Second valve (im parentheses) obdtamed from harmonization. 
Fic. 9. Values of (p+ q) obtained by breaking down the 


key points. Final values from harmonization are shown in 
parentheses for comparison. 
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be typical. A fairly 


other 
problems have been solved and the differences 
were approximately the same. There is, of course, 
no certainty that such will always be the case, 
and no claim is made that the errors will never 
be greater. But nearly the same situation exists 
in other approximate methods, such as the Ritz 


large number of 


method, which is widely used in problems of 
vibrations. For practical purposes the suggested 
procedure may be viewed as giving a good first 
approximation to the value of (p+gq) at an 
isolated point. 

The above application of the mean four-point 
influence value is, however, incidental. Its main 
value lies in providing a rational approach to the 
problem of determining a complete solution for 
the values of (p+gq) for the whole bar by the 
process of iteration or harmonization. To illus- 
trate this we again refer to Fig. 3 and observe 
that points 1, 2, 3, and 4 have been chosen in 
such a manner that they furnish key points for 
the direct transition to the finer mixed network 
shown in Fig. 9, this transition to be effected by 
means of the simple Liebman equation, 


(p+q).= > = dat Dot}. + >a), 


used either normally or diagonally at all points 
except those on the curved boundary. The final 
broken-down values of (p+q) are shown in Fig. 9. 


Improvement in the Transition from a 
Coarser to a Finer Network in a Mixed 
System 


At a curved boundary, where the values of the 
function change abruptly, it is desirable to have 
a much finer network than at straight portions 
where the changes in function are much smaller. 
To facilitate the work of iteration or harmoniza- 
tion Shortley and Weller! give the following 
improvement formulas: 


a= t(Watve tity. de 
Yr= i(Vitvaty, +y.), 
Y= &(2yotVst¥mt+2yi), 


¥i=9(2pat3pitynt3yx), 
Vi Ve (4b +5, +¥mtbnt5y;), 


the notation being as shown in Fig. 10. 
However, the introduction of a few auxiliary 
points (S), in the manner shown in Fig. 10, 
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Fic. 10. Sketch showing transition from coarse to finer 
network in a mixed system. 


eliminates the need for special equations and 
makes it possible to pass from one network to 
another twice or half the size of the first by the 
direct and simple use of the Liebman equation, 
employed either normally or diagonally. The 
values at the auxiliary points marked S can 
obviously be calculated by the Liebman equation 
used diagonally. Thus 


VS5s= (bo tYmtvnty;)- 


Points g’, b, g, 7, 7, k, etc., on the rows of the fine 
network adjacent to the coarse one can now be 
improved using the auxiliary points S as follows: 


Wo=i(Vsstysitys,+y.), 
Vo=4(Wor tVssty,t+v-), etc. 


The final values of (p+q) shown in Fig. 9 
have been obtained from the values at the key 
points by the above method. 


The Use of Conformal Transformations 


In problems of elasticity it often happens that 
the boundaries contain circular holes, grooves, or 
other discontinuities. In such problems, Car- 
tesian coordinates necessitate correction factors 
at the curved boundaries. Furthermore, since 
such discontinuities are generally sources of 
stress concentration, it is also necessary to use a 
much finer network near these sources, which 
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results in measurable increase in the amount of 
work. 

Considerable labor can be saved in such cases 
by combining the four-point influence method 
with polar coordinates. The use of such coor- 
dinates eliminates the need for correction factors, 
and gives from the start as large a number of 
points on the curved boundaries as conditions 
require. This automatically results in a finer 
network in the regions surrounding the sources of 
stress concentrations, and a progressively coarser 
network in the regions removed from the origin. 

The theoretical justification for this procedure 
rests on the possibility for a conformal trans- 
formation of rectangles into circles and vice 
versa, Fig. 11, by means of the equations 


u=logr (a); v=6 (bd). (4) 


It can be shown that, for analytic functions, 
all the essential relations between points in 
Cartesian coordinates, such as the nine-block 
procedure suggested by Shortley and Weller, are 
also valid at the corresponding points of the 
polar coordinate system. 












LEGEND: 


© — Known boundary values from 
photoelastic stress pattern, 








@ — Theoretical value; 





@ — Unknown interior points. 











Fic. 11. Drawing showing conformal transformation of a 
quadrant of a circle into squares. 
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Fic. 12. Photoelastic stress pattern of a ring subjected 
to diametral compression. Loads and dimensions as shown 
in Fig. 13. Model fringe value=500 psi tension or com- 
pression; thickness = 0.1795 in. 


Example 2.—<An ideal example illustrating the 
combined use of the four-point influence method 
and polar coordinates is the circular ring under 
diametral loads. Figure 12 shows the photo- 
elastic stress pattern of such a ring. The dimen- 
sions of the polar coordinates which are necessary 
in order to transform the ring into an approxi- 
mately square network are shown in Fig. 13. 
These dimensions Au and Av are calculated from 
Eqs. (4a and b). 

In order to have a square network it is obvi- 
ously necessary to divide the interval A—A, 
Fig. 13, into an integral number of equal inter- 
vals, say four. The value of Au representing one 
side of the transformed squares is given by 

Au= (log 2—log.1)/4=0.1735. 
But for a square network 
Av= Au. 
Hence 
Aé= Av=0.1735 = 10°, approximately. 

It is obvious that in cases such as rings it 
may not be possible to find a polar network 
which will transform into perfect squares. These 


errors can, however, be made small enough to be 
negligible. 
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Furthermore, in the immediate vicinity of the 
load, it is convenient to replace the concentrated 
load P by a radial stress system‘ 


o,= —2P cos 6/ art, (5) 


in which P is the applied load, r and @ the polar 
coordinates, and ¢ the thickness of the ring. 

Figure 14 shows a typical repeating section of 
the ring after transformation. It also shows the 
three key points as well as systems of coordinates 
in different directions for each key point, and the 
corresponding values of (p+ q) obtained from the 
stress pattern of Fig. 12. 

Figure 15 shows the three resulting influence 
curves and the mean influence values of (p+4q),, 
(p+q)2, and (p+q)3. These were broken down 
and transferred to the finer network shown in 
Fig. 16 by means of the Liebman equations used 
normally or diagonally. The results from the 
first breaking down of the key points are given 
by the upper and heavier number in each block. 
For purposes of comparison the values at con- 
vergence are shown in parenthesis in both Figs. 
9 and 16. 

The suggested procedure is therefore as follows: 
We determine several key values by the four- 
point influence method for a coarse network, 
break these down by use of the Liebman equation 





Fic. 13. Sketch showing key points in a quadrant of a ring. 


4 See S. Timoshenko, reference 3, p. 104. 
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to a finer network in the manner shown in Figs. 
9 and 16, and then test the solution by traversing. 
In the problems which we solved, one—or at 
most two—traverses were all that was needed. 
However, no general rule need be laid down for 
the number of traverses required as this is 
automatically indicated by the changes produced 
by the iteration. 


Discussion of Results 


Inspection of Fig. 9 shows some rather inter- 
esting results. The upper value in each block is 
obtained from the key points calculated from 
the mean four-point influence, whereas the lower 
values have been obtained from repeated travers- 
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Fic. 14. Drawing showing coordinates and their direc- 
tions for the key points in a repeating section of the con- 
formally transformed ring. 
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Point @: y, (mean ordinate) = 100 
100 from Harmonization 
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Fic. 15. Influence curves for the key points in the ring. 


ing or iteration and may be taken as the values 
at convergence. The closeness of the results is 
indeed striking. Thus, at point G5, the cor- 
responding values are 3060 and 3005, and at D5 
they are 3535 and 3500. 

Inspection of Fig. 16 shows an even greater 
degree of agreement between the unimproved 
values obtained from the key points and those at 
convergence obtained after many traverses. At 
the key points themselves these values are 
— 2520, —1255, and 100, compared with — 2420, 
—1280, and 100 at points (1), (2), and (3), 
respectively.® 

Limitation of space does not permit further 
illustrations. Other problems involving bipolar 
coordinates and the rational approach have, 
however, been worked out. In general, where the 
boundaries are curved, the mean unimproved 
influence values after a conformal transformation 
represent a closer approximation to the true 
values than those from the original boundaries. 

The original paper presented before the con- 
ference also contained a discussion of a funda- 


5 Contours or isopachic curves for the ring under con- 
sideration, based on the values from the numerical solution, 
and a composite photograph of the theoretical isopachic 
curves and those obtained by means of interference fringes 
from the surface of a steel model can be found in the paper 
given in footnote 3. 
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Fic. 16. Values of (p+q) obtained by breaking down the key point values. Final values 
from harmonization are shown in parentheses for comparison. 


This 
theorem and its applications have since been 
published, and will therefore be omitted from 
this paper.® 


theorem 


mental 


on stress trajectories. 


® See Max M. Frocht, Photoelasticity, Vol. I, pp. 205-213 
(John Wiley and Sons, Inc., New York, 1941). 
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Classification of Isotropic Points as Defined by 
o.-o:=0 within a Regular Region* 
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1. Limitations of the Proposed Problem 


S seen by the title of the present paper we 
are not going to consider isotropic points 
located on the boundary of the elastic body, nor 
shall we deal with singular points which cor- 
respond to applied external efforts, such as a 
force or a moment, concentrated at the point of 
application. However, the methods as developed 
in the present paper may be readily extended to 
cover the excluded cases. For the present, we 
assume that the stress components oz, oy, Tz, and 
the displacements u and v are one-valued finite 
regular functions of the coordinates x and y. 


2. One-Valued Determination of Principal 
Stresses 


Having two principal stresses and two letters 
a, and a2 to denote them by, let us establish a 
one-to-one correspondence by means of the 
condition 


O12 02. (2.1) 


If the principal stresses are different in value, the 
question as to which one is a; and which one is 
o2 is now answered in a unique sense by (2.1). If, 
however, the principal stresses are equal in 
value, the above question. becomes vacuous. 
Hence, (2.1) leads to a one-valued determination, 
of o; and a2 in any case. 

Considering the difference of the principal 
stresses ¢,— 02, we derive from (2.1) 


o1—022 0. (2.2) 
To compute that difference, we may use 
(o1—02)?=(¢,—0,)*+41:2," (2.3) 


in Cartesian coordinates, or 


* Paper presented at the Eleventh Semi-Annual Eastern 
Photoelasticity Conference, May 24-25, 1940, at the Car- 
negie Institute of Technology, Pittsburgh, Pennsylvania. 
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(o;— 02)? =(0,— 09)? +47,¢" (2.4) 


in polar coordinates. In either case, a considera- 
tion of (2.2) will lead to a one-valued deter- 
mination of the difference of the principal 
stresses. 


3. Isotropic Points 


Points at which the difference of principal 
stresses is equal to zero play an exceptional rdéle 
relative to stress trajectories, isoclinic and 
isochromatic curves, etc. All those curves show 
characteristic irregularities at such points, namely 
angular points, knots, multiple points, isolated 
points, etc. For this reason a special terminology 
is introduced and points at which 


o1—02=0 (3.1) 


are called isotropic. Considering the stresses, 
isotropic points do not display any irregular 
features, since the stresses are finite and can be 
represented by Taylor’s expansions. The in- 
fluence of isotropic points is structural: they 
determine the topological structure of the 
systems of stress trajectories, isoclinic and iso- 
chromatic curves, etc. 

An isotropic point, P, is called isolated if a 
circle of radius p>0O may be drawn around it 
such that P (the center) be the only isotropic 
point within that circle. 


4. Indefiniteness of ¢ at Isotropic Points 


The direction angle ¢ of a principal stress is 
given by 
tan 2@=2r,,/(¢,—a,) (4.1) 


in Cartesian coordinates, and by 
tan 2(@— 60) =27,¢/(0,—<9¢) (4.2) 


in polar coordinates. The above formulas become 
indefinite if, and only if, both numerator and 
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Fic. 1. 


denominator of either fraction vanish. In such a 
case (2.3) and (2.4) give 4, —o02=0, which means 
an isotropic point. Hence, at regular (non-iso- 
tropic) points the value of ¢ is definite and can 
be found from (4.1) or (4.2), while at isotropic 
points neither (4.1) nor (4.2) gives a definite 
answer for @ on account of the indeterminate 
form 0/0. 

The reason for the above failure of ¢ at iso- 
tropic points is that @ is discontinuous at such 
points, and thus no value at all is being assigned 
to ¢ at isotropic points. 


5. Index of an Isotropic Point 


Let P be an isolated isotropic point within a 
regular region. Let ABCDEFA be a circuit about 
P such that P is the only isotropic point within 
that circuit and no isotropic point is on the 
path of the circuit, either. Let us now consider 
the variation of ¢ as corresponding to a motion 
from A to A along the closed circuit ABCDEFA 
(Fig. 1). The variation of ¢ will be continuous 
since there is no isotropic point on the circuit. 
Back at A, the ultimate orientation of principal 
directions has to be congruent with the initial 
orientation at A. Moreover, the first and the 
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second principal directions, separately, have to 
arrive at A in positions congruent to their 
respective positions at start, because the prin- 
cipal directions cannot exchange their rdéles of 
being ‘‘first’’ and “‘second”’ during the circuit: 
such an exchange between o; and oe, which are 
continuous, would be possible only at a place 
where o;=¢2. But there is no place of this kind 
on the circuit—it would be an isotropic point, 
which we have excluded. ; 

Hence, the principal directions (01, o2) have 
rotated through an integer multiple of 180° during 
the circuit. This can be written as 


gue =jr, 


in which 7 is an integer (positive, negative, or 
zero, as depending on the intrinsic nature of the 
isotropic point P and independent of the shape 
of the circuit). The number 7 may be called the 
index of the isotropic point P: 


1 
j= pao. 
T 


For a regular (non-isotropic) point, gd@=0 
and thus j7=0. The converse is not true, since 
there exists an isotropic point (No. 8) for which 
j=0. As will be shown later, unit is an upper 
bound to j, hence all possible values of j are 
given by 


(5.1) 


(5.2) 


j=1,0, —1, —2, —3, «> (5.3) 


6. Airy’s Stress Function in the Vicinity of 
an Inner Point 


Airy’s stress function, F, is given by 


F=—r*+-—r’ cos 2(60— 69) 
2 2 





ymt2 | 
+> : cos (mé—X,,) 
™=lm+1bm 
Ma 
— cos ((m-+2)0— 1a) | (6.1) 
m+2 


In the above, T, S, 40, Lm, Mim, Xm, Mm are arbi- 
trary constants, and L,, 20, M,,2 0 for all values 
of m. 
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Absent from F are terms with negative or frac- 
tional exponents of r as well as the following 
terms: 


xlogr, ylogr, x6, y@, logr, 8, 
rlogr, 76, sin 26, cos 26, (6.2) 


because such terms would give rise to either 
infinite or many-valued displacements and 
stresses, which is excluded by the condition that 
the point P be an inner point of a regular region. 


7. Order of an Isotropic Point 


Using Airy’s stress function to compute the 
stresses in polar coordinates, we find that S=0 
for an isotropic point while the term with T has 
no influence whatsoever upon the formation of 
stress trajectories, isoclinic and isochromatic 
curves, etc., because the direction angle ¢ of the 
principal stresses does not depend on the term 
with JT. Hence, we may confine our consideration 
to that part of F which is given by the infinite 
series in (6.1). The value of m corresponding to 
the initial term in that series (first term actually 
present) is called the order of the isotropic point 
and is denoted by m. We may write formally, 


Order =n=minimum of m for 
(Limy Mm) #(0,0). (7.1) 


As the summation in the series (6.1) starts with 
m=1, the order, m, cannot be less than one, 
hence, the possible values are 


n=1,2, 3,4, -°:. (7.2) 


8. Isotropic Points without Complications 
An isotropic point of order 1 for which 
Li~M, (8.1) 


will be referred to as isotropic point without 
complications. If, to the contrary, L,=M,, the 
point will be called isotropic with complications. 
It will be shown that isotropic points without 
complications necessarily are tsolated points. 
Hence, isotropic lines are made up of isotropic 
points with complications. 


9. The Index-Order Equation 


For an isotropic point without complications, 
the initial term of Airy’s stress function deter- 
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mines completely and exhaustively the structure 
of stress trajectories, etc., at the isotropic point. 
To simplify the analysis, let us consider a one- 
term stress function and drop the superfluous 
subscript n: 


ynt2 i # 
F=— | cos (n@—d) 
n+itn 


M 
+—— cos ((n+2)0—1) (9.1) 


n+2 


in which 
L20, M20, L¥M. (9.2) 


Computing stresses from (9.1) and using them 
to determine ¢ by (4.2), we obtain by a subse- 
quent differentiation after rather lengthy iden- 
tical trigonometrical transformations, 


do 
2— = —(n—1) 
dé 
L?— M? 
+ --—— -— ——_———., (9.3) 
L?+M?+2LM cos (206—y"+A) 





Using (9.3) in combination with (5.2) to compute 
the index j, we obtain the index-order equation: 


j=—(n—1)+sgn (L—M), (9.4) 
in which Kronecker’s sign function, sgn x, is 


J 


1-—O 


| £34 5 


-2 eto 
~3 @e 








of 
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Fic. 2. 
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equal to one, zero, or minus one, according as its 
argument x is positive, zero or negative. Hence, 
considering (9.2), we note 


(1 if 
| —1 if 


L>M, 


sgn (L—M)= L<M 


(9.5) 
Figure 2 shows a graphical representation of the 
index-order equation, ring-shaped marks having 
been used for L>M, disks for L<M. The 
limitation n2 1 (7.2) now immediately gives rise 
to the limitation (5.3) of the index. 


10. The Isotropic Point of Index Zero 


Figure 2 shows the existence of such a point: 
j/=0, n=2, and the ring-shaped mark means 
L>wM. To have the simplest representative, let 
us put 


M=0, A=0, (10.1) 
which leads to 
, L 
F=—r' cos 20=~—(x'—y’), 
6 6 
¢,=0, o=2Lr*?cos26, 1,4=Lr* sin 28, 
(10.2) 
o,=—2Ly*, o,=2Lx?, 1,,=0, 
1=G,, O2=02, O1—02=2Lr*, o=0. 


The result ¢=0 reveals the most surprising thing, 
namely that the stress trajectories are rectilinear 
lines x=const. and lines y=const. The only 
interruption suffered is at (0,0) at the isotropic 
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Fic. 4. N=4. 


point. Except for this, there is no further struc- 
tural disturbance, and as shown in Fig. 3 we 
obtain a regular network of stress trajectories. 
In Fig. 3, as well as in all other figures which 
show stress trajectories (Figs. 4-6), we use a con- 
tinuous line for the first stress trajectories and a 
dotted (interrupted) pattern for the second stress 
trajectories. 

Isotropic points of the above kind have been 
actually observed on photoelastic photographs. 
Their fundamental ‘photographical’ charac- 
teristic is that (1), they are dark spots, cannot be 
removed by any rotation of the plane of polariza- 
tion of light, and (2), they do not produce a knot 
of the isoclinic system—an exception (the only 
one) from the general rule saying that an iso- 
tropic point is a knot of the isoclinic system. 
The stress trajectories are to be drawn straight 
through as if there were no isotropic point. 


Il. The Number R of Stress Rays 


In the vicinity of an isotropic point the system 
of stress trajectories is so shaped that some stress 
trajectories pass through that point. As stress 
trajectories usually have multiple points (angular 
points) at an isotropic point and appear broken 
and ramified, let us denote by R the number of 
rays (semi-lines) which are stress trajectories 
emanating from an isotropic point. For an ordi- 
nary, non-isotropic point the two orthogonal 
stress trajectories through it would give R=4 
rays. Omitting all of the mathematical analysis, 
we now state the resulting relative to R: the 
order n of an isotropic point being given, R may 
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Fic. 5. N=1 and N=12. 


have just one of the following three values: 
2n, or 2n+2, or 2n+4. (11.1) 


Particulars depend on the comparative sizes of 
L, M and n in (9.1), namely: 


(1), R=2n+4if L—M<O; 
(2), R=2nif (L—M)/(L+M) 
> 
1/(n+1); 6011.2) 
(3), any one of the three values (11.1) 
actually possible if 


0<(L—M)/(L+M) <1/(n+1). 





The Catalog Number N of an Isotropic 
Point 

For the sake of an orderly exposition of special 

isotropic points, let us introduce a _ catalog 

number N for each isotropic point as defined by 


N=}(11n+j—R-—2). (12.1) 


The above formula is arranged in such a way 
that the catalog numbers proceed in order of the 
natural sequence 


N=1, 2, 3,4, --- (12.2) 


without repetition and without omission. Any 
two isotropic points which differ either in n, or 
in j, or in R, will get two different catalog numbers 
by (12.1). The following table gives an exposition 
of the characteristics of isotropic points Nos. 
1-16: 
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( 12.3) 
13. Examples 


Isotropic points Nos. 1, 4, 5, 8, 12 and 16 are 
those most frequently observed in experimental 
photoelasticity. The figures as given below show 
the stress trajectories for those points. Some 
pairs of isotropic points, differing in order only, 
have identical stress trajectories. The difference 
in order would show up in the photograph of the 
isochromatics as the density of the fringe dis- 
tribution depends on n. Point No. 8 was shown 
in Fig. 3. It might be of interest to have a 














Fic. 6. N=5 and N=16. 


complete catalog up to, say, N=16. This work 
could be done by means of the formulas given in 
the present paper. For a more detailed account, 
containing some auxiliary formulas, see author's 
paper in Russian.! 


‘ Experimentalnie Methodi Opredelenia Naprajenii_i 
Deformazii v Uprugoi i Plasticheskoi Zonakh” (ONTI 
NKTP, Moscow, U.S. S. R., 1935), pp. 88-115. 
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Three-Dimensional Photoelasticity Using 
Scattered Light 


By R. WELLER 


Department of Mechanical Engineering, The State College of Washington, Pullman, Washington 


(Received February 8, 1941) 


It is possible to use the polarization caused by the scattering of light within a 
photoelastic model in place of the usual analyzer in a photoelastic polariscope. 
This enables one, in effect, to use an analyzer which lies within the interior of the 
model, and so make analyses of three-dimensional stress systems very conveniently. 
This 1s accomplished in practice by projecting into the model a beam of polarized 
light which has been collimated and passed through a slit, thus confining the beam 
to a thin sheet of light. Such a sheet will illuminate any desired plane section of 
the model, and in this section will appear interference fringes from which the 
stresses in the section may be determined. By successive examination of plane 
sections one may completely investigate the model. In this method the spacing of 
the fringes is the significant feature, rather than their order, as in previous methods. 
It is possible, in many cases, to determine the directions of the principal stresses 
at interior points as well as the maximum shear at such points. Boundary stresses 
are evaluated both as to direction and magnitude. Equipment for this method 
differs in constructional details from that employed in plane stress studies in that 
a source of light of higher intensity must be provided and a slit is used to restrict 


the beam. 





Introduction 


NTIL recently the photoelastic method for 

stress analysis has been limited to plane 
stress systems. There now exist, however, two 
truly three-dimensional methods which are 
theoretically capable of determining the stress 
system in any shape whatever. These are, in 
order of their discovery, the fixation method! and 


! The first mention of the possibility of this method 
occurs in J. C. Maxwell, ‘On the equilibrium of elastic 
solids,” Trans. Roy. Soc. Edinburgh 20, Part 1, 87 (1849— 
50). The model material used by Maxwell was apparently a 
gelatin obtained from the sturgeon. See also A. G. Solakian, 
“A new photoelastic material,’"” Mech. Eng. 57, 767, et seq. 
(1935). In German, see G. Oppel, ‘‘Polarisationoptische 
Untersuchung raumlicher Spannungs- und Dehnungzu- 
stande,”” Forsch. Geb. Ingenieur. 7, 240 et seg. (1936). 
(An English translation is available as Technical Memo- 
randum No. 824 of the National Advisory Committee for 
Aeronautics entitled ‘‘The photoelastic investigation of 
three-dimensional stress and strain conditions,” April, 
1937.) R. Hiltscher, ‘‘Polarizationoptische Untersuchung 
des raumlischen Spannungszustande im konvergenten 
Licht,” Forsch. Geb. Ingenieur. 9, 91, et seg. (1938) and 
A. Kuske, ‘“‘Das Kunstharz Phenolformaldehyd in der 
Spannungsoptik,”” Forsch. Geb. Ingenieur. 9, 139, et seq. 
(1938). In the United States, Hetenyi has discussed this 
method, e.g., ‘‘The fundamentals of three-dimensional 
photoelasticity,”’ J. App. Mech. 5, 149, et seg. (1938). 
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the scattered light method.? This paper considers 
the latter. 


The Stress-Optic Law 


The methods of photoelasticity are based on 
the fact that suitable transparent materials 
which possess optical isotropy when unstressed 
become birefringent under load. The develop- 
ment of this birefringence obeys the following 
relationships which are known together as the 
stress-optic law: 

1. Light passing through the stressed material 
divides into two components polarized at right 
angles to each other. These vibrate along the 
directions of the “effective principal stresses.’ * 
That is, if the stress ellipsoid at a. point is 
centrally intersected by the plane of the wave 


2 R. Weller, ‘‘A new method for photoelasticity in three 
dimensions,” J. App. Phys. 10, 266 (1939). 

3 —. C. Drucker and R. D. Mindlin (‘‘Stress analysis by 
three-dimensional photoelastic methods,’’ J. App. Phys. 11, 
724, et seq. (1940)) have used the term ‘‘secondary principal 
stresses.’’ Since secondary stresses have, ordinarily, a 
rather different meaning, it seems desirable to refer to 
them as “effective principal stresses.” 
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front, the section cut out will be an ellipse. The 
principal semi-axes of this ellipse are called the 
effective principal stresses in the plane of the 
wave front. 

2. The relative velocity of the two light com- 
ponents as they pass through such a point is 
proportional to the difference between the 
effective principal stresses. 

The relative velocity in (2) is a measure of the 
rate at which the phase shift between the two 
components increases along the light path. If we 
express this shift in wave-lengths (assuming 
monochromatic light) we may write: 


CdR/ds = S,—Sb, (1) 


where dR/ds is the rate of growth of the phase 
shift in wave-lengths per unit path length, S, 
and S, are the effective principal stresses, and 
C is a constant representing the stress-optical 
sensitivity of the model material. In the event 
that the direction of the original light beam 
happens to coincide with one of the three prin- 
cipal stresses at a point, then S, and S, will be 
the other two principal stresses; but in general 
this is not the case. 


Analysis by Scattering 


It is convenient to consider the scattering of 
light from the point of view of the absorption and 
re-emission of light by vibrating particles. The 
following are results of the simple theory: 

1. A beam of plane polarized light falling on a 
material particle in an isotropic medium will 
polarize it and cause it to vibrate in a direction 
parallel to the electric vector of the light, that is, 
in a direction normal to the “plane of polariza- 
tion”’ and lying in the plane of the wave front. 

2. The intensity of light scattered by a vi- 
brator in a given direction is proportional to the 
square of its apparent amplitude when viewed 
from that direction. By ‘‘apparent amplitude” 
is meant the component of amplitude normal to 
the line of sight. 

Now consider that a beam of initially plane 
polarized light passing through a stressed model 
is continuously analyzed by the. scattering 
process. It will be convenient to assume that the 
beam has been collimated and passed through a 
slit so that a plane section of the model is illu- 
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minated. At a given point within the piece, the 
induced motion of a scattering particle will be 
the result of the combined action of the two 
light components. This will generally result in 
the execution of an elliptical orbit in the plane 
of the wave front. Suppose, however, that the 
two components arrive exactly in phase. The 
ellipse then degenerates to a straight line at 45° 
to the two components and parallel to the direc- 
tion of the electric vector of the initial polariza- 
tion. If this point be observed from a point 
outside the model and along a line of sight 
parallel to this vibration, no scattered light will 
be observed. If, however, the two components 
arrive at this point exactly opposite in phase, 
the ellipse will degenerate into a straight line at 
right angles to the previous one and this will now 
be normal to the same line of sight. A maximum 
of intensity of scattered light will then be seen. 
Obviously if one should choose a direction of 
observation at right angles to the previous one, 
but still in the plane of the wave front, these 
effects will be reversed, the region which ex- 
hibited no scattering before will now provide a 
maximum and vice versa. 

Accordingly, as the light progresses through 
the loaded model and the two components 
progressively get in and out of step with each 
other there will appear variations in the intensity 
of scattered light along certain observation direc- 
tions. These will constitute interference fringes 
very similar to those observed in the conven- 
tional plane photoelastic effect. An exception to 
the above will be observed in the event that the 
observation takes place along the direction of 
vibration of one of the components, that is, along 
the direction of one of the effective principal 
stresses. In this case any scattered light will be 
the result of vibration due to but one of the 
light components, and quite naturally this cannot 
result in interference effects. Along any other 
direction the fringes will be observed, but with 
varying degrees of distinctness, the best direc- 
tions being those mentioned first. In fact the 
disappearance of the fringes when viewed along 
a direction of effective principal stress affords a 
convenient means for determining these stress 
directions. 

The spacing of the observed interference 
fringes will decrease with an increase in stress 
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intensity and a corresponding increase in the 
rate of growth of the phase shift. In Eq. (1) 
if the value of dR is taken as one wave-length, 
then the value of ds will be the path length along 
which one cycle of interference takes place; that 
is, it will be the spacing between two consecutive 


interference fringes in the direction of the 
original light beam. We may then write 
C/d=S,—Sh, (2) 


where d is the fringe spacing. This gives the 
average value of the effective stress difference 
between the two fringes. Unless the stress vari- 
ation is behaving badly this may be taken as the 
value at the midpoint. This is the more nearly 
true as the stress is increased and the fringe 
spacing decreases, the fringe spacing acting as 
the gage length over which the stress measure- 
ment is made. 

Drucker and Mindlin*® have pointed out that a 
correction must be made to Eq. (2) since a rota- 
tion of the axes of effective principal stress about 
the wave normal will cause a further crowding 





Fic. 1. Fringe pattern of notched tension model. Uniform 
tension in the shank is indicated by uniform fringe spacing. 
In the notch the fringes are more closely spaced at the 
boundaries indicating higher stress than in the center of 
the section. 


612 








Fic. 2. Microphotometric trace of notched section of 
tension model. In practice, the size of this trace is about 
ten to twenty times that of the negative. Measurements of 
fringe spacing may be easily made with a steel scale. 


of the fringes and lead to values of stress which 
are too high. If this rotation is small compared to 
the rate of relative retardation, the correction 
may be neglected in practice. 


Examples 


It will be convenient at this point to consider 
some examples of the above phenomena. The 
first of these is that of tension. Figure 1 shows 
the fringe pattern observed in a tension member 
of circular cross section having a semicircular 
notch. Here the direction of the initial light beam 
was normal to the side of the model and at right 
angles to the tension. The direction of polariza- 
tion made an angle of 45° to the direction of the 
tension in order that the light beam should be 
split into two equal components. Here Sa—Sp 
will be the actual tensile stress. The model was 
photographed from a direction lying in a plane 
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Fic. 3. Stress distribution across notched section of 
tension model. Values of spacing are obtained from Fig. 2 
and inserted in Eq. (2) to obtain stress values at mid- 
points between fringes. 





parallel to the wave front of the light and 
making an angle of 45° to the direction of the 
tension. This direction of observation is, of 
course, either parallel or perpendicular to the 
direction of the original polarization. In the 
event that it is parallel, the center of a dark 
fringe will lie at the point of entrance, and vice 
versa. This corresponds to the use of either 
crossed or parallel polarizers in the conventional 
plane-stress photoelastic procedure. 

Uniformity of tension in the shank of the 
model is evidenced by the fact that the inter- 
ference fringes are equally spaced, indicating a 
constant rate of phase shift between the two 
components of the light beam as they progress 
through the piece. The use of a uniform tension 
specimen affords a convenient method for cali- 
brating the model material. It is found, of 
course, that the value of C in Eq. (1) is simply 
the fringe value for a unit thickness in the con- 
ventional photoelastic method. For photoelastic 
Bakelite BT-61-893 this is approximately 86 lb. 
per sq.in. per fringe for direct stress, or half of 
this in terms of shear. 

The stress distribution across the notch is not 
uniform, as evidenced by the variations in the 
spacing of the fringes. In order to measure con- 
veniently the spacing of the fringes in such a case 
one may photograph the pattern and examine 
the negative with a microphotometer. Figure 2 
shows such a record of the fringe spacing in the 
notch of Fig. 2. Measurement of this micro- 
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photometric trace provides data for plotting the 
curve of Fig. 3, indicating the stress variation 
across the reduced section. 

A second example of the procedure is shown in 
Figs. 5 and 6. The model is a shaft subjected to 
torsion and having in it a Woodruff key slot. 
Figure 4 shows the pattern as viewed along the 
shaft axis, whereas Fig. 5 shows the view from 
the side. In each of these cases the illuminated 
section is normal to the direction of observation. 
The direction of polarization for Fig. 4 is at 
right angles to the shaft axis and in Fig. 5 it is 
parallel to the axis. This results in having the 
center of a light fringe lie in each case at the 
point of entrance of the light. The two sets of 
fringes are therefore actually different aspects 





Fic. 4. Fringe pattern due to torsion in shaft with key- 
way. This pattern is observed when the direction of obser- 
vation lies along the shaft axis. 





Fic. 5. Fringe pattern due to torsion in shaft with key- 
way. This is the shaft of Fig. 4 viewed from the side. The 
direction of polarization of the initial beam has also been 
rotated through ninety degrees. 
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of the same set. This is shown by Fig. 6, which 
gives the shape of the microphotometric traces 
for the two patterns. It will be observed that for 
all practical purposes the traces are identical. 

It is of interest to recall that the case of torsion 
was mentioned by Solakian in 1935' but received 
little attention since the existence of a fringe 
pattern in the chosen section was not easily 
explained. It should be pointed out that his 
pattern is actually correct if one considers the 
excessive deformation to which the model was 
subjected in this particular instance. It is ordi- 
narily considered that the principal stress direc- 
tions in such a case are the radii of the shaft and 
a set of 45° spiral helices winding about the shaft 
axis. In the case of large strains this is no longer 
true.‘ 

Another example of the utility of the scattering 
procedure lies in the fact that it may be applied 
to two-dimensional stress systems in such a way 
that the entire stress determination may be 
carried out optically. This has been treated by 
Drucker and Mindlin.’ 


Boundary Stresses 


In most applications of the photoelastic 
method one is interested chiefly in boundary 
stresses. Here the scattering method enables one 
to obtain complete information. At a free bound- 
ary one of the principal stresses lies normal to the 
boundary and is of zero magnitude. The other 


‘See B. R. Seth, ‘‘Finite strain in elastic problems,” 
Phil. Trans. Roy. Soc. A234, 231, et seg. (1935). Here it is 
shown that the angle departs from 45° by an amount 
6a/8 where @ is the angular strain and a the radius of the 
shaft. For @ of unity and a of 0.5 this gives 3.6° to be sub- 
tracted from 45°. The writer has been troubled by this 
phenomenon. If the direction of polarization of light enter- 
ing a circular shaft in which stresses are fixed, but which 
is examined by scattering, is parallel to a principal stress 
spiral at the point of entry, it remains so up to the centroid 
of the section. No fringes are therefore observed in the 
lower half of the shaft. As soon as the light passes the 
central plane, however, the direction of polarization is 
such as to give rise to fringes which are now parallel to 
the central plane and the effect is exactly as if the lower 
half of the model were cut away and light projected directly 
into the central plane. In fact, this cutting away was per- 
formed and the effect observed to be identical. The 
parallel fringes coincide with those observed from a direc- 
tion lateral to the shaft in the usual scattering procedure. 
Where the model is loaded in the usual way the strain is 
small and the effect mentioned above is hardly discernible. 
Traces of distortion may be observed in the upper half of 
torsion models, however, and these may never be com- 
pletely removed by any adjustment of the apparatus. 
But rotation of the polarizer to near coincidence with the 
stress directions will enhance the effect, even in cold 
loading. 
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FiG. 6. Microphotometric traces of torsion patterns. The 
right trace is obtained from Fig. 4, the left from Fig. 5. 
It will be observed that these are practically identical. 


two lie tangent to the boundary. If one projects 
the light beam into the model at right angles to 
the surface, the resulting fringe spacing will give 
the difference between the two tangent stresses. 
If the model is observed from various directions 
in a plane tangent to the boundary, there will be 
found two directions along which the fringes 
will vanish. These will be the directions of the 
tangent stresses. If the direction of one of these 
is known, the light may then be projected parallel 
to it. The fringe spacing under these conditions 
will give the difference between the other tangent 
stress and the (vanishing) normal stress, this 
difference being equal to the tangent stress 
itself. The second tangent principal stress may 
be determined similarly, or one may use the 
original difference obtained when the beam was 
projected normal to the boundary. 


The General Case 


If one desires to evaluate the stresses at an 
interior point he may proceed as follows. When 
light is projected through such a point at random 
the spacing of the fringes straddling the point 
will give the effective principal stress difference 
as defined previously. The minimum spacing will 
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Fic. 7(a). Schematic layout of polariscope for observation 
by scattered light. 


be obtained, however, when the light travels 
along the direction of the intermediate principal 
stress. The spacing will then be due to the dif- 
ference between the maximum and the minimum 
principal stress, that is, to the maximum shear. 
This may be obtained by a rotation of the model 
until the fringe spacing is a minimum, that is, by 
systematically searching for this position. When 
one has obtained the maximum shear and the 
direction of the intermediate principal stress in 
this way, the other two principal stress directions 
may be found by locating the two directions of 


observation along which the fringe system 
vanishes. Drucker and Mindlin* have suggested 
an alternative procedure in which one measures 
the fringe spacing from several directions and 
solves for the stress system by a consideration of 
these spacings and the angles of light projection. 
So far as the writer is aware neither of these 
suggestions has been put into practice. 


Apparatus 


The type of polariscope employed for plane 
stress investigations is not particularly suited to 
the scattering procedure. Equipment more satis- 
factory for this work is shown in Fig. 7. It is 
somewhat more compact than the two-dimen- 
sional apparatus and, in general, less expensive 
to construct. Several features of difference are 
as follows. The light beam is projected vertically 
upward in order to permit of observation always 
in a horizontal direction. The light source must 
be much more intense due to the inefficiency of 
the scattering process. No analyzing unit is 
required, but an immersion tank and fluid are 
necessary to eliminate boundary refraction. A 
slit must be provided to permit examination of 
single plane sections. 

The light source used has been the Type H-6 
high-intensity, water-cooled, 1000-watt mercury 
arc lamp, supplied by either the General Electric 
Company or the Westinghouse Company. Smaller 
lamps may be used but they make a dark room 
necessary and lengthen photographic exposure 
inconveniently. 














Fic. 7(b). Polariscope for observation by scattered light. The photograph on the left shows the instrument as constructed. 
The view at the right shows the parts before assembly. 
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Light from this lamp is collimated by a pair of 
condensing lenses, each 6.5 inches in diameter 
and having a focal length of 13 inches. The lamp 
is arranged so that it may be moved vertically 
and thus brought to the focal position of the 
lens combination. 

The light then passes through a Polaroid disk 
of diameter equal to that of the lenses. Provision 
is also made for the installation of a quarter- 
wave plate should this become desirable. Both 
Polaroid and quarter-wave plate may be rotated 
independently within the housing, each being 
provided with a graduated scale. 

The light then passes through a slit six inches 
long and of controllable width, the latter varying 
from about 3} inch to the closed position. The 
entire housing is arranged to rotate on its 
support, the slit always remaining parallel to the 
lamp burner for maximum light efficiency. The 
slit is of the two-jaw variety, both jaws moving 
simultaneously so that the center of the observed 
section does not alter with slit adjustment. 

The model to be observed is placed over the 
slit and light projected into it. Since the face of 
the model through which the light must enter 
is not, in general, plane and normal to the beam, 
the model must be immersed in a liquid having 
a refractive index equal to that of the model. 
One therefore requires a tank with transparent 
bottom and sides. Various liquids are possible in 
this connection. Best results have been obtained 
by the writer with “‘Halowax Immersion Fluid”’ 
which is water-white, of low viscosity, and com- 
paratively inexpensive. 

It does not seem possible to design a loading 
frame which will be suitable for any large variety 
of three-dimensional problems. Rather it has 
seemed advisable to build small loading fixtures 
for each investigation. Since the loading device 
must be partially or completely immersed, the 
use of fixtures of any size means a large tank and 
a large quantity of fluid. 

Photography may be accomplished without 
difficulty. During the writer’s first experiments 
he used a 100-watt mercury arc and found that 
exposures of several hours were necessary. The 
use of the 1000-watt lamp reduced the time to a 
few minutes. With a slit width of 0.1 inch, a lens 
having a speed of f/6.3, a film of Weston (day- 
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_ light) speed of 50, and a Corning No. 5 mercury 


green filter, the time required is about 5 minutes. 
Using a Wratten No. 77 filter the exposure is 
about 1 minute. Without the filters the exposure 
is approximately 10 seconds but the fringes are 
not so distinct. 


Conclusions 


Theoretically, there now exist two methods by 
which one may completely determine the stress 
system in any shape. The extent to which this 
will be true in practice depends on the develop- 
ment of the required laboratory techniques. It is 
the writer’s opinion that the accuracy obtainable 
from the scattered light method will be approxi- 
mately equal to that expected from the two- 
dimensional procedure. It is possible that the 
investigation of thin sections will require the use 
of a Babinet compensator, since there may be 
only one fringe or a fraction of a fringe visible. 

The scattering method seems to have the fol- 
lowing advantages over the fixation procedure: 
the gage length is reduced since, in most cases, 
the fringe spacing is less than the thickness of 
slices which might be cut out; the deformations 
are normal in amount, whereas in the fixation 
method they are about 25 times as great; the 
model is not destroyed and may be observed 
under progressive loading and unloading; para- 
sitic optical effects such as initial anisotropy may 
be observed easily. The last of these deserves 
special mention since it is not possible in three- 
dimensional work to subtract the stress due to 
initial conditions from those observed under load 
since the stress tensor may rotate during loading. 

On the debit side the scattering procedure 
requires a new form of apparatus for convenient 
work; the technique of making measurements is 
new and needs development; it is difficult, for 
example, to make measurements close to a 
boundary; the fringes are less distinct and 
become progressively more so as very thick 
sections are encountered. 

It is probable that the fixation and scattered 
light methods may, at times, be profitably used 
in combination, the model having the stresses 
“fixed” within it by heat treatment and being 
scanned with a scattering polariscope thereafter, 
or perhaps before being sliced. 
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A Photoelastic Study in Vibrations* 
BY W. M. MURRAY 


Massachusetts Institute of Technology, Cambridge, Massachusetts 


NDER the heading ‘‘Dynamic Stress Prob- 

lems’’ one includes all those cases in which 
the stresses vary with time as compared with 
static stresses which are independent of time, 
or almost so, if we neglect the effect of creep. 
In brief, the group of dynamic problems can be 
divided into two parts; first, those dealing with 
transient conditions of stress which die down in 
a short time, as in the case of impact, and 
secondly, those of a cyclical nature which re- 
peat themselves at regular intervals of time as 
in the cases of sustained vibration or repeated 
loading. 

To all these problems the photoelastic analysis 
is not applicable since, like other methods used 
for determining stresses, it possesses certain 
limitations. However, provided that these limita- 
tions are understood, the application of polarized 
light to the solution of dynamic problems is a 
powerful tool for the engineer. Until recently 
the photoelastic analysis has been applicable 
only to problems involving plane stress or to 
those three-dimensional problems which could 
be simplified to two dimensions, but means have 
now been developed for making photoelastic 
analyses of three-dimensional static stress prob- 
lems. It is very doubtful whether these methods 
can be applied to dynamic conditions, so that 
the study of this kind of problem is still limited 
to two dimensions. Even so, a wide field is open 
to the investigator. 

The development in photoelastic technique 
which has likely been of most importance in the 
study of dynamic problems was the introduction 
of highly sensitive materials (the synthetic 
resins) which facilitate the use of monochromatic 
light and photographic recording, thereby elimi- 
nating the tedious and uncertain observations 
made by matching colors with a compensator. 
Certain work! yielding satisfactory results has 


* Paper presented at the Eleventh Semi-Annual Eastern 
Photoelasticity Conference, May 24-25, 1940 at the Car- 
negie Institute of Technology, Pittsburgh, Pennsylvania. 
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been carried out on rotating disks with the use 
of white light, a celluloid model and a compen- 
sator, but the labor involved must have been 
very great even for this special study which was 
considerably reduced in complexity. This prob- 
lem of the rotating disk is on the border line 
between static and dynamic problems since it 
involves motion but no fluctuation in stress 
provided that the speed of rotation remains 
constant.? A similar problem but one which 
belongs truly in the realm of dynamic stresses is 
that of rotating gears which have been investi- 
gated.* This latter problem represents a case of 
cyclical or repeated stressing and as such lends 
itself to somewhat simpler type of analysis than 
the problems of transient stress. In this latter 
field notable work has been done in this country 
by Frocht,‘ and in Japan by Tuzi® and Sekiguti® 
who have investigated stress conditions set up 
by impact. 

Because of the rapidity with which these 
transient stress waves travel, the apparatus 
required to photograph and study them is 
somewhat complicated and for this reason the 
writer has confined his attention, for the present, 
to the study of cyclical stress fluctuation which 
can be examined with very much simpler 
equipment. Although the equipment and meth- 
ods used can be applied to a vast number of 
problems the writer has investigated in particular 
the stress condition set up in a cantilever beam 
under resonant vibration, since the available 
apparatus was readily suited to the problem. 

As is the case with most investigations, 
several arrangements of apparatus are usually 
tried out before the most satisfactory results 
can be obtained, and for this study the existing 
photoelasticity apparatus was adapted by re- 
moving the testing machine and replacing it 
with the vibrator and synchronized shutter. 
The model, a bar of transparent Bakelite with 
a rectangular cross section 0.4’ in depth and 
0.37”. wide, was clamped firmly (allowing some 
clearance) to a length of 3’ Duralumin I section 
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on the bottom of which was mounted an electro- 
magnet. By passing a 60-cycle alternating cur- 
rent through the coils of the magnet the whole 
system could be vibrated including the Bakelite 
bar, which, if its length were adjusted properly, 
would vibrate in resonance with an amplitude 


vr 


up to about }”. The 60-cycle vibration was 
much too rapid to be followed by eye, so a 
synchronized shutter consisting of a rotating 
paper disk with a slot operated by a synchronous 
motor taking power from the same source was 
introduced. This proved moderately 
satisfactory since by rotating the frame of the 
synchronous motor one was able to change the 
phase relation between the shutter and the 
vibrating beam and all positions of the cycle 
could be examined. The main disadvantages 
with this arrangement were its crudeness and 
the great length of time required to make a 
single picture and the fact that of three necessary 
variables, phase, frequency and amplitude, only 
the phase relationship of the shutter and the 


set-up 





Fic. 1. The electrical vibrator with Polaroid disks. 




















Fic. 2. Fringe pattern obtained with the electrical vibrator. 
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magnitude of the vibration could be controlled, 
the latter by regulating the power input to the 
electromagnet. 

It was recognized immediately that certain 
refinements would be necessary in order to make 
this apparatus more convenient with the result 
that a smaller length of the Duralumin I was 
selected, and upon this the magnet was mounted 
on the lower flange and on the upper flange a 
suitable clamp for the Bakelite beam was placed. 
To eliminate all the lenses of the previous 
photoelasticity apparatus the Duralumin beam 
with the electromagnet and Bakelite model was 
mounted on rubber pads in the wooden frame 
as shown in Fig. 1 and Polaroid disks were set 
up in front and behind the model. This arrange- 
ment was very simple and compact and could 
be adapted to any light source for either visual 
observation or photographic recording of the 
fringe pattern. As the light source a General 
Radio Strobotac or Strobolux was mounted 
behind the model. For visual observation this 
arrangement was excellent, since by suitable 
adjustment of the frequency of the flashing light 
source it was possible to make the photoelastic 
stress pattern appear to be motionless (Fig. 2) 
or to give the beam the appearance of very slow 
motion so that all positions of the cycle could 
be observed. For any visual examination the 
results were excellent. The apparatus, however, 
still lacked the desirable feature of a variable 
frequency which is necessary since the natural 
frequency of all the models cannot be conveni- 
ently arranged for 60 cycles per second. With 
the two requirements for more nearly mono- 
chromatic light and a variable frequency of 
vibration it was decided to go back to the 
regular apparatus which possessed a mono- 
chromatic light source and to develop a suitable 
vibrator and light source to go with it. 

For the vibrator, then, the necessary require- 
ments were that the frequency, phase and 
amplitude should all be variable and independent 
of one another. A miniature mechanical type of 
vibrator was chosen. This consists of two equal 
masses which can be rotated around parallel 
shafts in opposite directions and arranged in 
such a manner that their centers of gravity 
always occupy corresponding positions as indi- 


JOURNAL OF APPLIED PHYSICS 








aS" VT FT 








THEORY OF THE VIBRATOR 
THE HORIZONTAL COMPONENTS OF FORCE BALANCE OUT, THE VERTICAL 
COMPONENTS ADD UP. 


Fic. 3. 


cated in the slide (Fig. 3). Because of the fact 
that the mass and eccentricity for each weight 
are the same, the centrifugal forces will also be 
the same. With the correct relative orientation 
of the centers of gravity the components of 
force in one direction will balance each other 
out and the components in the direction at right 
angles will add up to produce a single force with 
a sinusoidal variation. This force is used to 
produce the vibration in the photoelastic model. 
Regulation of the fluctuating force must then 
be accomplished by controlling its magnitude 
and frequency. The former is adjusted by 
moving the centers of gravity of the weights 
nearer to or farther from their respective axes of 
rotation. Means have been provided whereby 
this may be done while the vibrator is in opera- 
tion so that one may increase or decrease the 
amplitude of vibration at will. In tests this 
feature has been found very useful and highly 
desirable. Figure 4 shows the vibrator set up in 
place with the model suspended below it; the 
rotating weights can be seen at the top but the 
variable speed motor which controls the fre- 
quency is not shown. The holder for the model 
is arranged so that it can be turned through 90 
degrees which facilitates use in different types 
of polariscopes. 

Originally a mechanical shutter was used in 
which the phase could be regulated relative to 
the model but this was not satisfactory because 
of the slack in the driving mechanism. Although 
some photographs were made with this equip- 
ment, it was discarded in favor of a flashing 
light source or mercury stroboscope of the 
Edgerton type. 
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Figure 5 shows a fringe photograph taken with 
this light source (and a Wratten No. 62 filter) 
which was mounted in the apparatus. The 
isochromatic lines are clear and well defined. 
In this case the model represents a cantilever 
beam 8” long with a depth of 0.322” and a 
breadth of 0.268’. The end held by the vibrator 
is enlarged to a depth of approximately 1}’’ and 
the change of section is accomplished with a 
fillet radius of 3’ to avoid any objectionable 
concentration of stress at this point. The 
frequency of the vibration was 3300 cycles per 
minute and the amplitude at the free end was 
0.56” although the photograph itself was taken 
when the displacement of the free end was only 
0.44’’. Although the result has been satisfactory 
in this case, one feels that it would have been 
better to photograph the conditions for maxi- 
mum displacement because at that instant the 
model is stopped instantaneously and there is 
less likelihood of blurring the picture because of 
the movement of the model. With a high speed 
light source this is not so important because 
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Fic. 4. Mechanical vibrator in the testing machine. 





Fic. 5. Fringe photograph obtained with the 
mechanical vibrator. 
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the duration of the flash may be just a few ten 
thousandths of a second. Compared with the 
electric the 
relatively permits longer time of 
exposure. 

Since a good deal of time has been given to 
the discussion of the apparatus and its develop- 
ment, let us now turn to the consideration of the 
results which can be obtained and how they 
correspond with theoretical computations. The 
following equation is given by Timoshenko,’ for 
the lateral displacement of a prismatical bar 
subjected to this type of vibration: 


stroboscope mechanical one is 


slow and 


i=w 


y=> X(A;i cos pit+B; sin pit) (1) 
i=1 


where y=the lateral displacement from the 
mean position. 

As the equation stands it represents the 
general expression for the displacement due to 
lateral vibration and takes account all 
possible normal vibrations. For the case at hand 
the bar is vibrating in the first mode only so 
that Eq. (1) may be written in the very much 
simpler form 


into 


y=X,(A1 cos pit+B, sin pit). (2) 


Hereafter, the subscripts 1 will be omitted as 
only the first mode is to be considered. The 
symbols on the right-hand side of Eq. (2) have 
the following significance : 





X =C,(cos kx+cosh kx)+C2(cos kx —cosh kx) 
+C;(sin kx+sinh kx) +C,(sin kx—sinh kx). 


A, B, Ci, Co, C3, and C, are constants to be 
found. 


p=a proportionality factor depending upon the 
frequency, 

= 2X frequency. 
k*= p*Dp/EIg. 
D=area of cross section of the bar. 
p=weight per unit volume of the material. 
E=modulus of elasticity. 
I=moment of inertia of the cross section. 
g= gravitational constant. 
L=length of the bar. 
x =distance measured along the bar. 
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Equation (2) contains seven constants which 
must now be evaluated according to the par- 
ticular conditions of the problem at hand. For 
the cantilever beam let x=0 at the fixed end 
and x=L at the free end. Then, one has the 
following conditions to fulfill: 


(a) Atx=0 y=0 hence (X)z.0=0 

(b) At x=0 slope=0 hence (dX /dx),.0=0 
(c) At x=L bending moment =0 hence (d?X /dx*),;.,=0 
(d) At x=L shear=0 hence (d°X /dx*),-,=0. 


From conditions (a) and (b) Cy=C3=0. 

From conditions (c) and (d) cos kL cosh kL 
= —1 which has for its first root RL =1.875 so 
that the numerical value of C, can be expressed 
in terms of C2 since 


cos kL+cosh kL 


C= —C- 
sin kL+sinh kL 





= —0.734 C.. 


Therefore 


y=C2{ (cos kx —cosh kx) 
—0.734(sin kx—sinh kx)} (3) 
X(A cos pt+B sin pt). 


If one measures the time ¢ from the mean 
position of the bar, then, when t=0, y=0, so 
that A=0 and one may write the equation for 
the displacement as 


y= BC2{ (cos kx —cosh kx) 
—0.734(sin kx—sinh kx)}sin pt. (4) 


The combined constant (BC2) of Eq. (4) can 
readily be found if one knows the displacement 
(other than 0) of some point along the beam and 
the corresponding value of pt. 

Equation (4) is now in form suitable for 
application to the specific problem of the model 
under test. However, before making any stress 
computations it will be well to find the effective 
length and to compare this with the actual 
length. 

From the previous equations 


k? = (1.875/Le)?= p(Dp/EIg)} 
so that 


3.515 ,Elg\ *} 
effective length = Le=| -(—) |. (5) 
p \Dp 
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The measured values for the beam are: 


L=length=8” 
h=depth =0.332” 
b=breadth = 0.268” 
frequency = 3300 cycles per min. 
E= 620,000 Ib./sq. in. 
p=0.050 Ib./cu. in. 
g=32.2 ft./sec./sec. 


Substitution of the numerical values in Eq. 
(5) gives for the effective length 


Le=8.21”. 


The measured amplitude of vibration at the 
free end of the beam was 0.56” (this corresponds 
to pt=7/2); therefore by substitution in Eq. (4) 


0.56 = (BC2) { [cos 1.875 —cosh 1.875 ] 
—0.734[sin 1.875 —sinh 1.875 ]}}, 


from which 
(BC2) = —0.280. 


y= —0.280} (cos kx —cosh kx) 
—0.734(sin kx—sinh kx)}sin pt. (6) 


When the fringe pattern was photographed, 
the displacement of the free end of the beam 
was 0.44” so that 


sin pt=0.44/0.56; 


8 
8 


therefore corresponding to the photograph 





y = —0.220{ (cos kx—cosh kx) 
—0.734(sin kx—sinh kx)}. (6a) 


Differentiating twice with respect to x gives 


M @y 
—— =—=+0.220 k?{ (cos kx+cosh kx) 
EI dx* 





—0.734(sin kx+sinh kx)} (7) 


where M=the bending moment and, since the 
stress at the extreme fiber is given by the 
equation 


f=+4Mh/2I, (8) 


f=+0.110k*Eh{ (cos kx+cosh kx) 
—0.734(sin kx+sinh kx)}. (9) 


The stress values along the edges of the 
beam have been computed from Eq. (9) and are 
shown plotted in Fig. 6. Two curves are given 
for the theoretical values; one is based on the 
actual length (8’’) of the beam and the other is 
based on the effective length (8.21). The third 
curve shows the photoelastic values obtained 
from the fringe patters (average of both sides) 
using a fringe constant at 87.7 lb. per in. per 
order. 
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The agreement between the actual length and 
the effective length of the beam is very good 
(within 3 percent). As one would expect, the 
effective length is slightly greater than actual 
length because of imperfect restraint at the fixed 
end and the fact that the vibrator was clamped 
on a little more than 8” from the free end. 

Since there is no question about the position 
of the free end of the beam, all quantities have 
been plotted with reference to it. The theoretical 
results are slightly higher than the photoelastic 
values but the agreement between the latter 
and the theoretical values computed using the 
effective length are very good indeed. It may be 
mentioned that the theoretical equations take 
no account of shear deflection nor of rotary 
inertia. However, since these quantities together 
will produce less than two percent difference in 
a beam of these proportions, it is reasonable to 
neglect them. 

In conclusion one may say that the photo- 
elastic analysis definitely substantiates the theory 
and establishes its validity under the test 
conditions. These are similar to those of certain 
types of magnetic fatigue testing machines as 


shown in Fig. 7. For anyone contemplating 

















Fic. 7. 
(Courtesy of A. V. de Forest, Massachusetts Institute 
of Technology.) 


research of a similar nature the author would 
suggest the use of some photographic method of 
recording the displacement of the beam and a 
high powered stroboscopic light source similar 
to that advocated by Dr. Sinclair.* 
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Abstracts 


HE following abstracts are from the Eleventh Eastern Photoelasticity 
Conference held at the Carnegie Institute of Technology, Pittsburgh, 
Pennsylvania, in May, 1940. Complete manuscripts for these papers were not 


available at the time of publication. 


A Photoelastic Investigation of the Stresses in the 
Foundation of the Sardis Dam Spillway. R. R. Puitippe,* 
U. S. Soils Laboratory, AND M. M. Frocut,} Carnegie 
Institute of Technology, Pittsburgh, Pennsylvania.—Intro- 
duction. Great progress has been made in recent years in 
the general field of stress analysis. Nevertheless, the 
structural engineer is still in search of a suitable method 
for the analysis of stresses in foundations. The theory of 
elasticity has yielded valuable results in a limited number 
of simple cases. However, where the boundaries become 
irregular and complicated, the theory becomes cumbersome 
and unwieldy. In these cases the photoelastic method of 
stress analysis has proved itself to be of special value to 
the structural engineer. Although some of the assumptions 
upon which the application of this method to soil mechanics 
rests may be open to debate, the practical value of the 
results as a first and useful approximation can hardly be 
doubted. The problem: A recent application of the photo- 
elastic method to the determination of boundary stresses 
in the foundation for the retaining walls in the Sardis 
Dam spillway is described. These walls are of the inverted, 
cantilever ““T’’ type which form the sides of the spillway 
and retain the overburden through which the spillway was 
excavated. Two typical sections of these walls were 
studied: one 41.0 ft. x30 ft. upstream of the control weir, 
and the other 48 ft. x44 ft. below the control weir. The 
purpose of this study was to evaluate the magnitude and 
distribution of normal and shearing stresses caused by 
combinations of overburden and hydrostatic loads judged 
to be most severe for these walls. An idea of the size and 
shape and loading of one of these walls can be obtained 
from Fig. 1. The method of analysis: Models were con- 
structed of the walls, the foundations of which were of 
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gelatin. The isoclinics and isochromatics were observed 
and recorded by sketches. Two methods were used for 
calculating the stresses: the “‘Filon Graphical Method’’!? 
along stress trajectories and the ‘‘Shear-Difference Meth- 
od’? along straight lines. The Filon graphical method 
became of questionable accuracy, due primarily to the 
difficulty of measuring the variables required by the 
method. On the other hand, the shear-difference method 
proved comparatively simple to handle and gave results 
which were well within engineering accuracy. The final 
distribution and magnitude of the normal and shear 
stresses are shown in Fig. 1. A discussion of the results is 
beyond the scope of this abstract. Gelatin models: The 
introduction of gelatin as a photoelastic material presents 
certain new problems in technique and calibration. The 
technique of the preparation of the model was discussed 
and three methods for the determination of the fringe 
value were given. The results from a uniformly distributed 
load over a finite region on the edge of a semi-infinite 
region were compared with those from a concentrated 
load on the same plate as well as with the linearly dis- 
tributed body compressions in a rectangular block of 
uniform cross section. Attention was called to the necessity 
for rigid temperature control of the gelatin model. The 
model for the Sardis Dam foundation was prepared and 
kept at about 65°F throughout the investigation. The 
effect of the body forces upon the stress pattern was also 
discussed. 

* Engineer, U.S. Soils Laboratory, Pittsburgh, Pennsylvania. 

+ Associate Professor of Mechanics, Carnegie Institute of Technology, 
Pittsburgh, Pennsylvania. 

1See Max M. Frocht, Photoelasticity (John Wiley and Sons, Inc., 
New York, 1941) Vol. I, Chapter 9. 

2See E. G. Coker and L. N. G. Filon, Treatise on Photoelasticily 
(Cambridge Press, 1930). 


3 See Chapter 8 of reference 1. _ 
4 See Figs. 10.16 and 10.17 of reference 1. 


The Annealing of Photoelastic Bakelite and the Re- 
sulting Effect on Its Physical and Optical Properties. 
G. H. LEE* anp M. Go.pstTEIn, Cornell University, Ithaca, 
New York.—In recent tests conducted in the Cornell 
University laboratory on annealed photoelastic Bakelite, 
certain inconsistencies in the properties of the material 
were noted. In an effort to determine the cause of these 
variations, an investigation was instituted to determine 
the effect of annealing on the properties of Bakelite, i.e., 
physical and optical creep, Poisson’s ratio, Young’s 
modulus and the stress and strain optical coefficients. 
The study of the annealing methods was made with the 
aid of a specially constructed annealing oven and a 
straining frame with a constant temperature testing 
chamber. The former was so designed as to permit photo- 
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elastic observation of a Bakelite sample by means of 
appropriate window construction in the top and bottom of 
the oven, as its annealing progressed. Care was taken to 
test all samples at the same temperature in order to 
eliminate any effect that may arise due to variation in 
temperature. The sample was fitted with a longitudinal 
and a lateral extensometer, and the photoelastic effect was 
observed with the aid of a Soleil-Babinet compensator. 
The procedure was to choose several unannealed, ground 
and polished plates of Bakelite (BT-—61-—893) and to cut 
them into a number of samples of equal size. Groups 
were then formed, each consisting of one sample from 
each plate. The samples from each group were then tested 
after having been subjected to the same annealing pro- 
cedure. The results of the given annealing were then taken 
as the average of the results of the tests on each group. 
Several groups having been formed, each was subjected 
to a different annealing procedure and the results com- 
pared. One group was tested in the unannealed condition 
for purposes of comparison with annealed specimens. The 
results of the study are shown by the accompanying 
tables. Note that two series of tests were made; one to 
determine the effect on properties if the samples were left 
in the oven for a time only sufficient to remove completely 
the initial stress, and a second to study the effect on 
properties due to varied time of annealing. It is desirable 
that the material used in photoelastic investigation have 
the following properties:—(a) low stress optical coefficient, 
(b) low creep (physical and optical). With this in mind, 
if Tables I and II be studied, it will be seen that the best 


TABLE I. Time of annealing just sufficient to remove initial stress. 


| | 





Younc's STRESS- STRAIN- | | 
ANNEALING | MopuULus Optica | Opticat | Unit | Optica 
TEMPERA- OF COEFFI- CoerFFI- | CREEP | CREEP 
TURE ELASTICITY) CIENT | CIENT (23 HR.) | (23 HR.) 
(°F) (LB./IN.?) |(LB./IN.FR.)| (FR.IN./IN.)| (IN./IN.) (FR.) 
Unannealed | 618,000 87.8 1980 0.00083 1.11 
184 632,000 88.8 | 2040 0.00104 | 1.41 
225 612,000 86.2 } 2050 0.00105 | 1.49 
242 602,000 | 85.3 | 2010 0.00103 1.49 
601,000 85.9 | 2020 0.00143 1.90 
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TABLE II. Effect of time of annealing on the properties of Bakelite 


temperature of oven 242°F. 
Younc's STRESS- STRAIN- 
Moputus | OPTICAL OPTICAL Unit | OPTICAL 
ANNEALING OF COEFFI- COEFFI- CREEP | CREEP 
TIME ELASTICITY CIENT CIENT | (23 HR.) | (23 HR.) 
CHR.) (LB./IN.?) |(LB./IN.FR.)| (FR.IN./IN.)| (IN./IN.) (FR.) 
0.44 602,000 | 85.3 | 2010 0.00103 1.49 
5.0 | 602,000 86.3 2088 0.00131 | 1.88 
12.0 | 618,000 89.1 | 2070 0.00146 1.80 
' 


compromise in these properties is obtained for the case of 
annealing temperatures between 240°F and 250°F, and 
for a length of time just sufficient completely to remove 
all initial stress. Accordingly, as the final result of this 
study on the effect of annealing procedure on photoelastic 
Bakelite, the following annealing procedure (for specimens 
of usual thickness) is recommended by the authors:— 
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(1) Preheat oven to 240-250°F, (2) place sample in oven 
for one hour, (3) oven cool the sample at a rate of ap- 
proximately 15°F per hour. 


* Now at Carnegie Institute of Technology, Pittsburgh, Pennsylvania. 


Some Examples of Failure Due to Stress Concentration. 
R. E. PETERSON, Westinghouse Electric and Manufacturing 
Company, East Pittsburgh, Pennsylvania.—The results of 
photoelastic investigations are used in industry to a greater 
extent than is generally realized. Stress concentration 
factors usually appear in engineering department design 
manuals years before becoming available in books on 
machine design. A considerable number of machine 
manufacturing companies have photoelastic equipment, 
ranging from simple set-ups to complete laboratories. 
The reason for all of this is that failures due to actual 
breaking of parts can in some cases be very dangerous and 
expensive, and therefore a considerable expenditure can 
be justified to improve design proportions. The chairman 
of the committee for this meeting asked that some examples 
of failure be described so that those engaged in photo- 
elasticity at universities might better appreciate potential 
applications of their work. In this abstract only a few 
typical examples will be given. Figure 1 shows fatigue 
cracks at the base of gear teeth, located exactly where 
photoelastic tests would show that the cracks should occur. 
The paper by Professor Dolan deals with this question 
in more detail. Figure 2 shows fatigue cracks started at 
the end of a keyway in an overhung street-car motor drive 
shaft. The variable stress was primarily bending. A similar 
case where the variable stress was primarily torsion is 
shown in Fig. 3. This failure was due to a resonance 
condition occurring in a large motor-generator set. Some 
work has been done on the general keyway problem using 
three-dimensional photoelasticity. The general field of 
fastenings should prove fruitful to the photoelastician. 
In Fig. 4 is shown a failure through the cheek of a large 
Diesel engine crankshaft due to variable bending. A 
similar case due to torsion is shown in Fig. 5. The problem 
of optimum crankshaft design is one in which photo- 
elasticity can be of considerable assistance, although, of 
course, other important factors must be kept in mind. 
Figure 6 shows a failure of the flanged shaft of a Diesel 
driven generator. A number of helical cracks (usually 
associated with variable torsion) have joined to produce 
the saw tooth appearance. If the fillet were larger in 
proportion to the shaft diameter, one large helical crack 
would form. If the fillet is very small compared to the 
shaft diameter (as in Fig. 5) the crack apparently is 
confined to the fillet proper. The rational design of filleted 
members calls for considerable aid from photoelasticity. 
Some twenty examples similar to the above were described 
in the presentation. Speaking for the industrial laboratories, 
we are very much interested in the photoelastic work being 
done at the universities and are particularly hopeful that 
the newer techniques will be given major attention. 


1T. J. Dolan, J. App. Phys. 12, 584 (1941) (this issue). 
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Frictional Phenomena. II 
A. Gases 


By ANDREW GEMANT 


The Detroit Edison Company, Detroit, Michigan 


Chapter II. Theory and Experimental Facts 
on Gaseous Viscosity 


N the basis of the kinetic theory, an explana- 

tion is given both of the viscous forces 
within a gas and of those external frictional 
forces operative between a gas and a wall. The 
extent to which the theory is able to reproduce 
the existing experimental data on the viscosity 
of different gases, and the variation of viscosity 
with pressure and temperature are shown. 


8. Kinetic Theory 


Turning our attention first to the viscosity of 
gases, it is instructive to review the experimental 
facts together with the corresponding theoretical 
considerations. The theory of gaseous viscosity 
is comparatively simple, and was developed 
along with the experimental research, whereas 
with liquids the theory still lags behind the 
experiments. The reason for this is the relatively 
simple structure of gases which led to the 
establishment of the kinetic theory. 

The kinetic theory allows certain equations to 
be deduced connecting the viscosity constant 7 
(Section 4) with other characteristic quantities 
of the gas. These equations were tested on the 
experimental material available. The agreement 
was found to be only partially satisfactory, and 
this led, in the course of developments, to con- 
tinual improvements of the theoretical formulae. 
A short review of this situation is now presented. 

We refer to Eq. (6), according to which the 
frictional force X, acting on the x-z plane in the 
direction of the x axis is given by 


X,= —ndu/dy, 


if the flow u depends only on y (see Fig. 3(a)). 
This type of flow is called a parallel flow, since 
all velocities are parallel to one direction (the 
Xx axis). 
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In gases molecules are so far apart from each 
other that the motion of the molecules is of 
primary. importance, and their mutual at- 
traction, which becomes essential in liquids and 
solids, is only of secondary, though not negligible, 
interest. The mechanism of the frictional forces, 
too, is based on the movement of the molecules. 

The mechanical explanation of the viscous 
force will be understandable if we recall one of 
the elementary laws of mechanics. The mo- 
mentum, M, of a mass my with velocity v is 
defined as myv. If a force f acts upon the body, 
we have the relation 


f=modv/dt 
or 


f=dM /dt. (14) 


In other words: The force acting upon a body is 
given by the increase of its momentum per 
second. 

This law can be used directly in this case. 
The viscous drag acting upon the x-z plane is 
due to the gas molecules impinging upon that 
surface, and trying, so to speak, to carry it 
along with them. The magnitude of the force is, 
according to (14), equal to the momentum 
transported through that surface in the unit of 
time. The amount of this momentum can easily 
be computed. 

Let the number of molecules per unit volume 
be m and the mass of each molecule be m. As to 
the velocity factor in M, the thermal components 
can be disregarded, since they are equal at both 
sides of the surface, and a finite transport of 
momentum is not possible. The speed component 
that has to enter in M is the flow velocity 


u=uytydu/dy, 


and here again only the second term is of im- 
portance. If the y component of the thermal 
motion of individual molecules is now denoted 
by v, and the number of molecules in the volume 
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unit having such a velocity component is de- 
noted by n,, then the number of such molecules 
hitting the surface unit of the x-z plane in the 
time unit is n,v, and we have 


dM du 
—=F nol my), (15) 
dt v dy 


summing over all the different v values. Here y is 
determined by the location of the molecule at its 
last collision before passing the plane x-z, since 
it had acquired the velocity u during that 
collision. Suppose a molecule hitting the x-z 
plane at O (Fig. 5) has had its last collision at 
O’, with a path r and a y component —y between 
the two points. We then have 


v/c=—y/r, 


where c is its thermal velocity. In reality c is 
distributed according to Maxwell’s law, but we 
operate here with its average value. Placing y 
into (15), we have 


1 du 
f=-—-— > n.mv"r. 


cdy v 


The average of 7 is called mean free path in 
the kinetic theory of gases, and is denoted by X. 
Thus, we write 

\ du 


f=-—-— E n,m’. 
cdy » 


The sum }>°,.”,mv* is shown in the theory to be 








Fic. 5. Diagram of molecule transferring momentum. 
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the pressure p of the gas, thus 
f=—(Ap/c)du/dy. 


Comparing the last equation with Eq. (6), we 
have 
n=Ap/c. (16) 


Using the relation, as supplied by the theory, 
for p: 


p=nmc/3, 
we also have for 7 the alternative equation: 
n=nmdc/3=0.33 prc, (16a) 


where p=nm=density of the gas. 

The physical reason for » being proportional 
to pdc is this: Both the density (p) and the 
temperature (c) favor an increase in the number 
of molecules which effect a transportation of 
momentum. Because of a finite gradient du/dy, 
the transported momentum per molecule is larger 
the farther away from the surface the last 
collision has occurred, that is, the larger the 
mean free path. 

Equations (16) and (16a) refer the viscosity 
constant 7 to other known molecular constants 
of the gas. It is a theoretical relation of remark- 
able simplicity which allows a check in the light 
of experimental facts. 

The factor } in (16a) is only an approximation, 
since a constant velocity was used instead of 
taking account of a Maxwell distribution. Several 
later calculations, especially by Chapman! have 
shown that the factor 0.33 has to be increased. 
Accordingly, we write 


n=0.49pdc. (16b) 


Equation (16b) can be used to compute values 
of the mean free path \ from measurements of 
the viscosity. The density of the gas is known, 
and the average velocity c is computed by means 
of the relation 

mc? = (8/mr)kT, (17) 


(k= Boltzmann constant = 1.37 X10~"* c.g.s. and 
T =absolute temperature). 

After \ has been computed, it is possible to 
make use of a further relationship to calculate 
the diameter o of the molecules, a quantity 
that enters the kinetic theory when we picture 
the molecules as rigid elastic spheres from the 
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standpoint of collision. The theory then gives 


A=1/v24r0n0", (18) 


which allows o to be calculated from \, and thus 
from viscosity data. The calculated values of \ 
and o are in fair agreement with values of the 
mean free path and molecular diameter as com- 
puted from entirely different experimental data. 

The following table gives the results of some 


computations, compiled by Kennard,’ based on. 


the data of several authors, chiefly Trautz*® and 
co-workers. The c.g.s. unit of the viscosity is 
called 1 poise. 


TABLE |. Mean free path and molecular diameter of gases 
computed from viscosity data (p=1 atmos., temperature 


= 15°C). 

Gas 7 (10°* Potse) | A (10-6 cM)! @ (1075 cM) 
Ammonia 97 4.5 4.4 
Carbon dioxide 145 4.2 4.6 
Ethylene 100 3.6 4.9 
Helium 194 18.6 a4 
Hydrogen 87 11.8 y 
Mercury vapor 219°C 470 8.3 4.3 
Neon 309 13.2 2.6 
Nitrogen 173 6.3 c ® 
Oxygen 200 6.8 3.6 
Water vapor 93 4.2 4.6 
Xenon 224 3.8 4.8 


9. Experimental Facts 


We now turn to certain fundamental experi- 
mental facts and compare them with the above 
theory. 

Substituting (18) into (16b) it can be seen 
that the density cancels out, since m is propor- 
tional to p. It follows that the viscosity is inde- 
pendent of the density, as well as of the pressure, 
if the temperature is constant. This independence 
is a consequence of » being proportional to the 
product (pA). It seems strange, perhaps, but it 
is confirmed, at least approximately, by the 
experiment. The viscosity of gases increases only 
slowly at high densities. The increase of the 
viscosity of air at 30°C up to 320 atmospherics‘ 
is shown in Fig. 6. The viscosity of gases appears 
to decrease at very low densities; a discussion of 
this effect will be given below. 

This independence of the pressure has also a 
practical importance. In certain cases—for in- 
stance, in the measurement of internal damping 
of solids by means of vibrations—it is desired to 
climinate the losses due to the viscosity of the air. 
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Fic. 6. Viscosity of air at 30°C up to 320 atmospherics 
according to Moulton-Beuschlein. 


The specimen is therefore placed in a container 
with reduced pressure. It has to be borne in mind 
that the losses due to a viscous flow of the air 
will not be reduced by such procedure. In case 
there is a turbulent flow (a nonstationary process 
as compared with the viscous flow), losses might 
be reduced if the pressure decreases. At high 
frequencies losses due to sound radiation will be 
reduced also; but this has nothing to do with 
the viscosity of the medium. 

The effect of temperature upon viscosity can 
be seen from (16b) in combination with (17). 
Accordingly, » ought to be proportional to \/T. 
Figure 7 shows, according to Kennard, the vis- 
cosity of some gases as a function of the absolute 
temperature. It can be seen that the increase is 
decidedly stronger than according to the 3 power, 
especially at higher temperatures where the 
slopes are nearly constant. 

Sutherland has developed a theory that ac- 
counts for this deviation. He introduces attrac- 
tion forces with a potential U(r) between the 
molecules. The effect of such forces is illustrated 
in Fig. 8. Molecule A, if no such forces are 
present, will pass by molecule B without collision 
(path 1). If the potential U is finite, A will 
either be deflected (path 2), or even driven 
against molecule B, causing a collision (path 3). 
For a given U, it will depend on the velocity of 
the molecule whether it will follow a path like 2 
or 3, slower velocities favoring stronger de- 
flections. 

In carrying out this idea quantitatively, 
Sutherland arrives at an effective cross section 
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F1G.-7. Viscosity of different gases as a function of 
temperature. 


Set, given by an equation: 
Serr = So(1—4U(e) /mce’*), (19) 


indicating an increase as compared with So, the 
real cross section, since U is negative. Substi- 
tuting Eq. (19) into (18) and (16), it is easy to 
see that the viscosity will take the form 


n=al'/(1+C/T). (20) 


This formula indicates an increase with tem- 
perature stronger than with the 3 power. The 
reason for this behavior is an increase of the mean 
free path with increasing temperature (path 2 
preferred to path 3). 

As a quantitative check indicates, even the 
Sutherland formula fails to account completely 
for the (7) curves. The Sutherland constant C 
in Eq. (20) was determined by Ribaud and 
Vasilesco® for air, No, A and COs. It increases 
slightly in the temperature range 0 to 1600°C. 
Some recent papers by Massey and Mohr® 
introduce a quantum-mechanical concept for the 
collision with fairly good success in computing 
viscosity values. Hellund and Uehling’ also de- 
veloped a general theory of transport in gases, 
using quantum-mechanical conceptions of the 
collision process. The equations for viscosity 
show the same dependence on temperature and 
pressure as the classical formulae, except that 
the coefficients contain the quantum-mechanical 
data. 

Itterbeek and Paemel*® measured the viscosity 
of argon at low temperatures, but could not 
reproduce the results by means of a theoretical 
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equation. Keesom and Keesom® find a linear 
increase of viscosity of hydrogen with tempera- 
ture in the range of 20 absolute. 

As to the viscosity of chemically different 
gases, Eqs. (16) and (18) indicate a propor- 
tionality to m'/o*. Since the cross section in- 
creases with mass, a large variation of 7 will not 
be expected. Still, the variation of m is larger 
than that of ¢, and, therefore, the viscosity is 
expected to increase with the molecular weight. 
Both these conclusions are confirmed by the 
data as shown in Table I. 

Itterbeek and Claes” found the viscosity ratio 
deuterium /hydrogen very near v2, in agreement 
with the theory. With decreasing temperature 
the ratio decreased. 





_ ~| 
A 


B 2 
3 


Fic. 8. Collision of molecules under forces of attraction. 


Much research has been done on the viscosity 
of gaseous mixtures. Adzumi'! measured the 
viscosities of mixtures of low molecular weight 
hydrocarbons between 20 and 100°C, and deter- 
mined the constants in Eq. (20). Frequently 
maxima and minima in the curves referring the 
viscosity to the concentration ratio of the 
mixtures were observed. Schréer™ criticizes the 
work of Adzumi whose results do not agree with 
the theoretical formula of Schroer. 

The viscosity, as well as the thermal con- 
ductivity, decreases if a magnetic field is applied. 
Results of Senftleben and Gladisch'* refer to 
oxygen and nitrogen monoxide, and can be for- 
mulated by a few simple relations. 


10. Theory of Slip 


Up to this point, theory and experimental data 
on the viscosity constant 7 of gases were dis- 
cussed, without reference to the technique of 
measurements. Although this latter aspect will 
not be treated in detail, the principle of such 
tests has to be mentioned for a better discussion 
of slip. More details on measurements will be 
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given in connection with liquids (see Chapter IV) 
and solids. 

One of the most frequently used methods is 
that utilizing the flow in narrow circular tubes. 
With a given flow the pressure difference at both 
ends is measured, and these data allow the 
constant 7 to be computed. A simple apparatus 
of moderate accuracy on this basis has been 
recently described by Weber.'* Another typical 
method is that of two concentric cylinders, the 
outer being rotated at a constant speed, the 
inner elastically suspended. The torque of the 
inner cylinder, corresponding to a given speed of 
the outer, is a measure of the viscosity of the 
gas between the cylinders. Bearden'® made use 
of this method for a precision determination of 
the viscosity of air at 20.00°C. It is 1819.20 
+0.06 10-7 poise, calculated after Inglis'® had 
computed the effect of a minute eccentricity of 
the two cylinders. 

In both cases there will be a certain interaction 
between the gas and the solid walls, and this 
interaction is likely to affect the data. If such 
interaction could be avoided entirely, then 7 
could be measured without interference on the 
part of the walls. Such method is supplied by 
the determination of damping of freely propa- 
gating acoustic waves (see Chapter III). The two 
principles mentioned above are, however, more 
practicable, and, therefore, generally used. For 
this reason it will be necessary to consider the 
interaction between gases and walls; in other 
words, to consider the external viscosity of gases, 
as compared with the internal viscosity ». We 
shall see that the practical consequences of such 
considerations become important, especially at 
low gas pressures. 

Whereas 7, the internal viscous force was 
defined, according to Eq. (6), as proportional to 
the velocity gradient, an external viscous force 
between a gas and a wall must evidently be pro- 
portional to the velocity difference between the 
two, and, in case the wall is stationary, simply 
to the velocity u of the gas. The viscous force 
per surface unit acting between the gas and the 
solid is then given by eu, where e is the constant 
of external viscosity. 

Let us consider the layer of the gas nearest 
the wall. There will be two opposite forces acting 
upon it: The solid wall pulls it backwards by a 
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force eu, and the next gaseous layer pulls it 
forward with an intensity ndu/dz (z=perpen- 
dicular to the wall). Thus in equilibrium 


eu = ndu/dsz. (21) 


Introducing the so-called coefficient of slip ¢ by 
means of the definition 


f=n/e, (22) 
we have 
u=(du/dz. 


(23) 


Equation (23) is then to be used as the boundary 
condition in any mathematical treatment of a 
gaseous flow. 

An expression for ¢ (or e) on the basis of the 
kinetic theory can be given in the following 
way.” Again considering the layer next to the 
wall, it transfers momentum from the bulk of 
the gas to the wall, and the two evidently must 
be equal. As shown above (Eq. (14)), the viscous 
drag ndu/dz is equal to the transfer of momentum 
per second. Since this transfer in the bulk occurs 
in both directions, and in equal amounts, the 
transfer in one direction from the bulk to the 
boundary layer is 

Sndu/ dz. 


Now as to the transfer from the layer to the wall. 
The theory shows that } (mc) molecules strike the 
surface unit per second. The total momentum 
striking the wall is }ucmu. The specific assump- 
tion should now be added that the total mo- 
mentum is absorbed by the wall, none being 
reflected. This happens if the wall is rough, and 
the molecules are severally reflected to and fro 
before their final escape, thus losing their velocity 
component « completely. Equalling the two 
momentums, we have 


3ndu/dz=incmu 
and, comparing with Eq. (23): 
¢=2n/ncm. 


Substituting Eq. (16b) into the expression for ¢, 
and writing 1.0 instead of 0.98, we have 
t=. (24) 


Since \ is inversely proportional to the pressure, 
(24) can be given the form 


6 =f0/pP, (24a) 


JOURNAL OF APPLIED PHYSICS 








2) 


3) 





where ¢o is a constant independent of the pres- 
sure. It can be seen that the slip coefficient 
increases with decreasing pressure. 

It can be shown that its value at atmospheric 
pressure is too low to be of much importance. 
Consider Eq. (23) and write orders of magnitude, 
instead of actual figures. For du/dz, then, we 
write u,/Z4 where u,~velocity of the gas and 
Z,~dimension of the apparatus (radius of tube, 
for instance). We have, then, with Eq. (24) and 
Uw = velocity of gas at the wall: 


Uw /Ug~d/ Za. (25) 


\ at atmospheric pressure is ~10~> cm; hence 
Uw/U,gK1, in other words, slip is negligible. The 
boundary condition is u,,=0 with a fair degree of 
approximation. 

With decreasing pressure the mean free path 
increases, and the slip velocity cannot be neg- 
lected. This will result in a relatively increased 
flow, the same as if the viscosity constant were 
to decrease. Finally, when \>2,, the character 
of the flow changes entirely and is controlled by 
the external friction e, which is proportional to 
the pressure. This is the key for the understand- 
ing of the so-called “molecular air pump,” as 
shown in Section 11. 

For an illustration of the significance of a 
finite slip let us consider the flow, discussed in 
Section 7, which has been treated with the 
assumption that there is no slip at the boundaries. 
Let us now generalize the case and use Eq. (23) 
for a boundary condition. Whereas (11) remains 
unchanged, we have for u — 


u=(p/2n)(h?—2?)+¢ph/n. (11a’) 


Thus we see that to the parabolic component of u 
is being added another term, constant over the 
whole cross section. For the total flow V we 
have: 


V =2ph®/3n+2¢ph?/n. (11b’) 


In utilizing Eq. (12), we get for the total heat 
developed W, 


W =2p*h®/3n+2¢p*h?/n. (11d’) 
However, if we compute the heat developed 
within the liquid, using (11c) as before, we 


arrive at (11d) which is only the first term in 
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(11d’). This seeming contradiction is explained 
by interpreting the second term in (11d’) as the 
heat developed in the external friction process. 
That this conclusion is correct can be shown 
directly. Applying our general equation A ;= — uf, 
from Section 2, the viscous drag at the wall is 
—euy; thus, the heat developed is eu?,,, for both 
walls together 2eu?,,. Substituting u. from Eq. 
(11a’): 
Uw =ph/n, 


we obtain exactly the second term in Eq. (11d’). 


11. Experiments on Slip 


Before discussing some experimental facts, re- 
lated to the slip coefficient, a brief derivation of 
the flow formula for circular tubes should be 
given. The formula also will be used extensively 
throughout this whole monograph. 





Fic. 9. Cross section of cylindrical tube. 


Figure 9 is a cross section of a cylindrical tube 
of radius a and length /, at the ends of which the 
gas pressure is p; and pe, respectively. A steady 
stream of gas will flow through the cylinder, with 
a volume output of V per second. The gas 
velocity u in axial direction will be a function 
of the radius 7, as shown below. Considering a 


. cylindrical volume of radius r, thickness dx 


(x=axis), and applying Eq. (8) for the x com- 
ponent of the forces, we have 


ar-dp+2xrdxndu/dr=0, 
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and, integrating over the whole length: 


du r Pi-pe2 


dr 2n 1 


Integrating over r with the boundary condition 
(23), we have 


1 pi- 
u=—(a*—r*?+2fa) 


4n 


(26) 


In carrying out the integration 2x J/j’rudr over 
the cross section, the volume V is obtained: 


ma‘ 4°\ pi-— po 
yet ( gS yor 
8n a l 


With ¢=0, this equation reduces to the well- 
known Poiseuille formula. In the second term ¢ 
has to be replaced by {o/p according to Eq. (24c), 
in order to have V as an explicit function of the 
pressure. Eq. (27) is only valid if (pi—pe) is 
small as compared with p; otherwise in the 
integration of the second term we should have 
taken into account that ¢ is a function of the 
pressure. In the second term of the exact equa- 
tion the expression In (p;/p2) should stand in- 
stead of (pi— po) /p. 

Eq. (27) clearly shows, as pointed out in the 
previous section, that for normal pressures, when 
\/a<1, slip can be neglected. In reducing the 
pressure, the output increases, even if Ap = pi— pe 
remains constant. For very low pressures the flow 
becomes inversely proportional to the pressure 
(for given Ap). 

These predictions of the theory were verified 
experimentally by Stacy'® and van Dyke.'® It 
should be recalled that in deducing (24) the 
absorption coefficient of the wall (f) has been 
taken as unity. If the value of f is left open, one 
has for ¢ instead of (24) 


g¢=(2-f)a/f. 


. Millikan has calculated f from the measurements 
just mentioned, utilizing the rotating cylinder 
apparatus and using an equation in principle 
analogous to (27). The value of f was in certain 
cases 1.0, in others around 0.9, thus justifying 
the use of Eq. (24). 

Adzumi”® found that the slip coefficient for a 
gas in metal capillaries varies with the nature of 


(27) 


(24b) 
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the metal. He considers the gas molecules as 
being adsorbed at the wall, and derives ex- 
pressions for the slip coefficient on this basis. 
Important research in this field has been done 
by Knudsen and Gaede.” Let us consider Eq. 
(27) and see how the flow depends on the pres- 
sure, if the pressure difference is a given fraction 
of the total pressure, or (p1— p2)/p=constant. 
The first term is then proportional to p, the 
second becomes independent of p. If not the 
volume, but the total amount of the gas (volume 
times pressure) is asked for, then the first term 
becomes proportional to p?, the second to p. 
Thus, in plotting pV as a function of p (Fig. 10), 
the quadratic curve 1 would be obtained if 
Poiseuille’s equation should hold for the whole 
pressure range. But, since for lower pressures the 
function must be linear, curve 2 should be 
expected. In actual experiments curve 2 was 


pV 





Pp 





Fic. 10. Output of flow as a function of pressure in the 
Poiseuille and free-molecular range. 


verified in the lowest range, the transition, 
however, occurring, via curve 3, for reasons not 
quite clarified. 

The lowest range is of particular interest. It 
occurs when \A>a; thus, the first term in (27) 
can be neglected entirely. This type of flow is 
sometimes called the free-molecular flow, since 
intermolecular collisions do not occur at all. 
The only interaction is that between individual 
molecules and the wall. Gaede has utilized the 
flow in this range in his so-called molecular pump 
using rotating cylinders. One can see from curve 1 


of Fig. 10 that for a given fraction of pressure 
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difference (say 10 percent of the existing pres- 
sure), the output of a pump would be extremely 
low on the basis of Poiseuille’s law. But, on the 
basis of 3, which actually takes place, the output 
of a pump can still be appreciable, even when 
the absolute value of the pressure drops to a 
very low value. This is just the mechanism of 
the molecular pump. One cylinder rotates against 
another and, when the pressure has dropped to a 
low value, carries off the remainder of the gas by 


means of the external friction, operative in 
this range. 

An interesting application of slip at atmos- 
pheric pressure is Carson’s determination™ of the 
pore size in papers by using Eq. (27). He argues 
that-a in papers is so small that the second term 
cannot be neglected. In carrying out flow meas- 
urements across paper at 1 and 0.9 atmospheric 
the average pore radius is computed. The values 
vary from 0.2 to 1.2. 
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Velocity Measurement of Transient Mechanical Motions 
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Measurements of velocity of rapid non-recurrent mechanical motions by means of an 
oscillograph are described. Two methods of measurement are given: One by modifying the 
optical system of the oscillograph, and the other by employing a slotted disk mounted on 
mechanism under test. These methods offer practical means of time-distance relation measure- 


ments of small devices, such as spring-actuated levers, cam mechanisms, triggers, etc., whose 
performance would be affected due to the additional load required to operate auxiliary meas- 


uring instruments. 


HERE are many transient non-recurrent 

mechanical motions, whose time-distance 
relationship (velocity, acceleration) is difficult to 
measure unless elaborate apparatus is employed. 
Velocity of such devices as spring-actuated 
levers, cam mechanisms, triggers, over-center 
mechanisms and other “‘snap’’ mechanical move- 
ments, which occur in a very short interval of 
time and whose traveled distance is relatively 
small, may be easily measured by means of an 
oscillograph, which may be found in many 
electrical laboratories. 

An oscillograph is usually associated with the 
measurement of electrical phenomena, such as 
current, voltage, power and their phase relation- 
ship. Or it is widely used in the investigations of 
minute mechanical motions producing changes in 
the sensitively-tuned electrical circuits, such as 
measurement or recording of sound or of me- 
chanical vibrations in terms of electrical oscil- 
lations, by means of a microphone. 

Mechanical movements may also be studied 
directly, without having their motion relayed 
through the medium of an electrical device 
(contacts, microphone, photoelectric cell, etc.). 
In many cases this is of paramount importance, 
since the coupling of the mechanism under study 
to a device that operates the oscillograph vibrator 
involves absorption of some energy otherwise 
available for the operation of the device proper, 
or it creates the possibility of unaccountable 
delay in the time element recording, with conse- 
quent erroneous results. 

Take, for example, a small current-interrupting 
mechanism, whose contacts, when tripped, open 
with a quick motion, essential for successful arc 


634 


rupturing. The velocity and the acceleration of 
contacts separation are very important elements 
in the determination of the device’s performance. 

To measure the average velocity of the device’s 
travel is a relatively easy task, for the time (7) 
required to complete the travel could be meas- 
ured with two oscillograph galvanometers, or by 
other simple means, the traveled distance (.S) 
may be easily calculated, and the average ve- 
locity (V) may then be computed from the basic 
formula 


V=S/T. (1) 


However, it is well known that a mechanism, 
when started from rest, cannot reach its maxi- 
mum velocity instantaneously due to inertia of 
parts, and in order to determine the true charac- 
ter of the motion, it is necessary to know the 
instantaneous values of velocity during any 
interval of time of movement’s travel. 

If a point is moving along any path, the 
instantaneous velocity (v) is the average velocity 
over indefinitely small period of time, including 
the instant, or, 


v=ds /dt, (2) 


where ds is the distance traveled during the 
indefinitely small interval of time dt. Conse- 
quently, if the instantaneous positions of the 
mechanism could be recorded simultaneously 
with the elapsed time, the rate change in travel 
may be calculated and represented graphically. 
The problem, therefore, is to record the relative 
positions of the mechanism from start to end of 
its travel, either at regular and known distance 
intervals, or continuously, in relation to a known 
time scale. Thus, if a trace of the mechanical 
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motion could be produced simultaneously with a 
record of a timing wave (alternating current of 
known frequency) the instantaneous values of 
velocity may be calculated. 

If the oscillatory mechanical motion is small, 
of an order of 1 to 5 degrees, an oscillograph 
vibrator mirror may be fastened to the shaft of 
the moving element. The mirror, being very small 
will not materially affect the inertia of the 
mechanism under consideration. The mechanism 
is then mounted within the oscillograph case in 
place of a vibrator, and so adjusted that the 
movement of the mirror causes the reflected light 
beam to trace a record on the moving oscillograph 
film, thus recording the mechanism’s position at 
any instant of its travel. 

In case the mechanical motion is relatively 
large (about 60°) then the above method could 
not be used, inasmuch as the light beam would be 
thrown beyond the width of the oscillograph 
film. To diminish the scale of the beam throw the 
optical system of the oscillograph has to be 
modified by the use of a prism, located in the 
light beam path. 

1, it will be observed that 
the prism, which is attached to the shaft of the 
oscillating mechanism under investigation, is 
placed in the path of the light beam reflected 
from a galvanometer mirror. This galvanometer 
is not connected electrically, and thus the beam 
its mirror reflects is stationary. When the prism 
oscillates, the refracted beam produces a straight- 


By referring to Fig. 


PRISM ATTACHED TO 


MECHANISM UNDER 
INVESTIGATION 


Fic. 1. Schematic diagram of oscillo- 
graph elements showing the position of 
prism attached to the mechanism under 
investigation for prism method. 


Pee. 














Fic. 2. Oscillograph and auxiliary apparatus showing the 
prism attached to the mechanism and located in the path 
of the reflected light from the stationary galvanometer 
mirror (prism method). 


line trace on the oscillograph film, which when 
set in motion at right angles to the beam travel, 
will produce a curve in right-angle coordinates. 
Actual apparatus set-up is shown in Fig. 2. A 
small prism is mounted on a shaft rotating in fixed 
bearings. This assembly is placed within the 
oscillograph case so that the prism lies in the 
light beam path reflected from a stationary 
galvanometer mirror. The shaft carrying the 
prism is attached to the shaft of the oscillating 
mechanism under investigation. The mechanism 
may be mounted outside of the oscillograph light- 
tight enclosure. Since the prism is rigidly, con- 
nected to the mechanism, its movement must 
necessarily correspond to that of the mechanism. 
Consequently, the beam deflection, as recorded 
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on the film, will be proportional to the movement 
of the mechanism. 

To determine the position of the mechanism 
during any instant, the prism travel has to be 
calibrated. This is accomplished with a pointer 
attached to the prism shaft and a graduated dial. 
Curve plotted with beam deflection as abscissae 
and angular travel of prism as ordinates, makes it 
possible to determine the instantaneous angular 
positions of the mechanism from the trace on the 
film. 

An alternating current wave of known fre- 
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Fic. 3. Oscillogram showing a record 
produced by motion of the prism 
which was attached to the moving 
element (yoke of an ‘‘over-center”’ 
mechanism). Simultaneous record of 
time scale is made to permit calcula- 
tion of mechanism velocity (prism 
method). 


quency recorded simultaneously with the phe- 
nomenon on the same film furnishes the timing 
scale for determination of the elapsed time 
between various points of mechanism travel. 
From these data, the values of ds and dt may be 
obtained and substituted in Eq. (2) to calculate 
velocity (v) at any instant of the mechanism 
travel. 

An oscillogram of such measurement is repro- 
duced on Fig. 3. It represents a curve produced 
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Fic. 5. Curve showing relation between the angular 
travel and the velocity of the moving element (yoke of 
an “over-center’’ mechanism). 
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Fic. 6. Oscillogram showing the 
record of the angular travel of the 
moving element (yoke) of an ‘‘over- 
center’ movement. The record of an 
interrupted beam was produced by a 
slotted disk which was placed in the 
path of a light beam reflected from 
a stationary galvanometer mirror. 
Current (0.15 amp.) was carried by 
contacts attached to the yoke 
(slotted disk method). 


by an “over-center’ mechanism, in which a 
spring, when loaded, releases a yoke which 
swings through an angle of 60°. A space-time 
curve, Fig. 4 shows graphically the time rate of 
mechanism motion. From the data of this curve 
the relation between velocity and mechanism 
travel was determined and plotted, on Fig. 5. 
Velocity of mechanisms which do not lend 
themselves readily to practical means of fastening 





Fic. 7. 
mechanism under investigation, carrying a slotted disk, 
and mounted within the oscillograph case. The motor in 
foreground operated the mechanism through gear reducing 
means (slotted disk method). 


Apparatus showing the oscillograph and the 
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of prisms to their moving elements may be 
determined by a slotted-disk method. A light 
weight disk, with slots pierced radially and 
spaced at known angular intervals, is mounted on 
the spring-driven intermittently moving shaft of 
the mechanism under investigation. The mecha- 
nism is mounted within the oscillograph case so 
that the slots permit the light reflected from a 
stationary galvanometer mirror to fall upon the 
moving film during the disk’s rotation. This 
produces an interrupted trace on the film. As in 
the previous method, the timing wave is simul- 
taneously recorded on the film, to provide the 
time scale for calculation of the velocity. An 
oscillogram obtained by this method is repro- 
duced on Fig. 6. It is quite apparent that since 
the angular travel and the time elapsed may be 
measured, the velocity may be readily calculated 
as described above. The apparatus arrangement 
is shown on Fig. 7. 

It must be borne in mind that for measure- 
ments described above it is desirable to have an 
oscillograph with accessories permitting recording 
at relatively high film speeds (up to 35 feet per 
second).! If the transient is of very short duration 
(0.001 second) such equipment is essential. 

It has been found that these methods of 
measuring the velocity of transient mechanical 
movements offer results of relatively high degree 
of accuracy with the apparatus usually found in 
electrical laboratory, and are applicable to 
small mechanisms, whose performance would be 
affected by the additional load required to oper- 
ate the auxiliary velocity recording or measuring 
instruments. 

1S. Keilien, ‘“Oscillograph—versatile tool in studying 
transient phenomena,” Elec. J. 35, 61 (1938); S. Keilien, 


“Improving products through research and test for 
quality,”’ Product Engineering 10, 24 (1939). 
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A Simple and Rigorous Method for the Determination of the Heat Requirements of 
Simple Intermittently Heated Exterior Walls 


ELMER G. SMITH 
Agricultural and Mechanical College of Texas, College Station, Texas 


(Received December 17, 1940) 


Buildings such as churches and auditoriums that are not continuously heated need larger 
heating plants than continuously heated buildings because of the large transient heat currents 
that must be supplied to the cold walls while the temperature of the building is being raised. 
The present unsatisfactory practice is to estimate the magnitude of these currents by various 
unreliable empirical methods. The purpose of this paper is twofold. First, for simple exterior 
walls, an equation is derived by means of which tables can be computed that will make rapid 
and reliable determinations possible. Second, the equation is put into such a form that the 
actual calculation of the tables can be carried out by N. Y. A. boys whose services can be 
obtained without cost. The unidirectional heat flow equation is solved for the boundary 
conditions that the outdoor temperature remains constant, the temperature of the wall is 
constant previous to time ¢=0, and after that time a constant heat current is applied to the 
wall. The solution is obtained in a form simple enough for the boys to handle by putting part 
of it in dimensionless form which makes it possible to plot a family of curves, which in turn 
make it possible to obtain numerical solutions without resorting to the very tedious ‘‘cut 





” 
and try” process. 


N a building such as a church or an auditorium 

that is heated occasionally and allowed. to 
cool when not in use, the required capacity of the 
heating plant is determined, in part, by the mag- 
nitude of the transient heat currents into the 
walls and furniture during the time in which 
the temperature of the interior of the building 
is being raised. To date, no method for deter- 
mining the magnitude of these heat currents has 
been devised that is both simple enough and 
accurate enough to permit the rapid and reliable 
calculations that those engaged in this particular 
industry need to be able to make. It is the 
purpose of this paper to derive a rigorous equa- 
tion for a simple exterior wall and put it in a 
usable form. 

Since the heating plant is usually run at a 
rather high and constant rate while the tem- 
perature of the building is being raised, the 
assumption that heat is being delivered to the 
walls at a constant rate during the transient 
period is probably reasonably accurate. 

Choose the origin and signs as indicated in 
Fig. 1. Choose the temperature scale arbitrarily 
in such a manner that its zero corresponds with 
the initial temperature of the building and its 
contents. Let 


6=temperature on the arbitrary scale. 
x =distance from the indoor surface of the wall. 
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L=the thickness of the wall. 

i=the heat current per unit area at any point 
in the wall. 

r=the thermal resistivity of the wall. 

c=the thermal capacitance per unit volume of 
the wall. 

{=time. 


The differential equations for the wall then are 


— (00/dx)dx = ridx 
— (01/dx)dx =cdx(d0/ dt). 


The solution of the above equations is 


0=A, exp [—(rcp)'x ]+A.2 exp [+(rcp)'x ] 
i=(cp/r)'A, exp [ —(rcp)'x ] 
—(cp/r)'Ae exp [+ (rcp) ‘x ], 














+0 
A 
Interior of on of 
Building he oors 
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> ~_ +X 
L 
Fic. 1. 
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where A; and A» are independent of x and 
p=(0/dt). 
The boundary conditions are: 


i=1,1 when x=0 
0=R,i, when x=L, 
where 


i,=the constant heat current per unit area that 
is assumed to enter the indoor surface of 
the wall continuously after ‘=0. 

1 = Heaviside’s unit function. 


-R, [exp [+(rcp)'(L—x) ]+exp [- ware Gre 


8=——_ eee eee ee eee 


R, =the surface resistance per unit area of the 
outdoor surface of the wall. 

i, =the heat current per unit area through the 
outdoor surface of the wall. 


It should be noticed that although 7, is a con- 
stant, 7, is a function of time and must be 
eliminated before A; and A» can be determined. 

After A, and A» have been determined, the 
equation for the temperature at any point in 
the wall may be written 


Arexp [+(rcp)'(L—x) ]—exp C— (rep)(~L—x) 1) 





RoA(cp/r)Cexp [ + (rcp)L]—exp [= (rep)!LY]+Lexp [+ (rep)'L]+exp [— (rep)'LT] 


Application of the Bromwich integral to the above equation gives 








Ry | ex A)ML— > —(rcr)'(L—x 
ts tial =)J-texp te Loexp ex) —x) ]—exp [—(rcd)'(L—x) ]] 
@=— $- a ae ik Pitan alls 0 APB ' 
2aj (Ry (cr r)iCex xp [+(rcd)!L ]—exp [—(rex)'LJ]+[exp [+(rcd)!L |+exp C— (ran IL] : 


where j=(—1)!. To get rid of the radicals let A 


=—(u?/rc). 


When the above substitution has been made and account taken of the fact that \ goes around its 
contour twice when u goes around its contour once, the result, after simplification and reduction to 


trigonometric form, is 


9=— 





2rj 


In Eq. (1) and throughout the paper the broken 
line bracket () encloses an angle. It is evident 
that a double pole occurs at 
u=0 

and simple poles occur where 

Ryu sin (Lu)—r cos (Lu) =0. (2) 
The residue at the double pole can easily be 
found to be 

[Re +r(L—<x)]/r. (3) 

Before the final results of this investigation can 
be expressed as a formula, as they must be if 
they are to be useful, some easy way to solve 
Eq. (2) must be found. The method adopted can, 
however, be limited in range to those values of 
R,, L, and r, that occur in the walls of buildings. 


For the present, it will be assumed that such a 
method has been found. Let 


¢=any solution of Eq. (2). 
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rts rs exp (—u*t/rc)_Ryu cos ((L—x)u)+r sin (L—=)u)] 
ete else é ~de 
u*(— R,-u sin (Lu)+r cos (Lu)) 


(1) 


The Laurent e expansion about a: any of the points $ 
can now be found by the following process. Find 
the Taylor’s series about the point @ for each of 
the following quantities. 


1/u?, (4) 
exp (—w*t/rc), (5) 
Ryu cos ((L—x)u)+r sin ((L—x)u), (6) 
— R,-u sin (Lu)+r cos (Lu). (7) 


Multiply the expansions of (4), (5), and (6), 
together. In the expansion of (7) delete every 
term that contains 


R,-@ sin (Lo) —r cos (L¢) 


as a factor because by the definition of ¢ as a 
solution of Eq. (2) all such terms are zero for 
any possible value of ¢. Finally, divide the 
product of the expansions of (4), (5), and (6) by 
the remaining terms of the expansion of (7). 
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When this has been done the residue for any of the simple poles can be seen to be 


exp (—¢*t/rc)[ (Reo cos (¢(L—x))+r sin (@(L—x)) ] 


(8) 
¢*| — R, (L¢ cos (Lo) +sin (L¢)) —rL sin (L¢) ] 
Since the solution is the sum of the residues (3) and (8) multiplied by 7i,, the final result is 
| +2 exp (—¢*t/rc) Rw cos (¢(L—x))+r sin (¢(L—x)) | 
6=1,, Ry tr(L—x)—-r>d - ; (9) 
| > —«*R,-(L¢ cos (Lo) +sin (L¢)) +rL sin (L¢)) 
where ¢ takes on all the real values given by the following equation : 

Ry @ sin (Lo) —r cos (Lo) = 0. (10) 


The quantity to the right of the summation sign in Eq. (9) can be simplified to some extent. If the 
quantities cos (¢(1—-x)) and sin (¢(L—x)) are expanded in terms of the sines and cosines of (L@) 
and (gx), it will be noticed that the coefficients of sin (¢x) add up to zero because of Eq. (10). 
Also, from Eq. (10) cos (L¢@)=(R,-/r)¢ sin (Ld). When this substitution has been made and the 
result simplified, Eq. (9) may be written 


+2 exp (—¢"t/rc)(Rid +r) cos (gx) 
0=1,, Ry +r(L—x)—-r > . ) cos L 
| —_ o°(LR* ¢’*+rR _+Lr’) 


(11) 


By setting x=0 in Eq. (11), the expression for the temperature of the interior surface of the wall is 
obtained. If, in addition, the arbitrary temperature scale is replaced by the ordinary Fahrenheit or 





centigrade scale, the result may be written 


0,—0;=14 Re t+rL—r Dv 


| , a 


where 


0;= the initial temperature of the building and its 
contents. 

0, =the instantaneous temperature of the interior 
surface of the wall. 


In Eq. (12), the lower limit of the summation 
has been changed from —~« to 0, and com- 
pensated by introducing a factor of 2. This is 
possible because ¢ occurs only as ¢ and according 
to kq. (10), for every possible positive value of 
@ there is a negative value of equal absolute 
magnitude and vice versa. In order to put Eq. 
(12) in terms of the instantaneous temperature 
of the air inside the building, apply Ohm’s law 
to the interior surface of the wall, which gives 
Reis =Oa—9s, (13) 
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+2 CX 


¢?(LR? ¢*+rR +Lr*) |’ 


p (—¢*/rc)2(Rve +r)| 


where 

R,=the surface resistance per unit area of the 
indoor surface of the wall. 

0,=the instantaneous temperature of the air 
inside the building. 


Substitution of the value of 6, from Eq. (13) into 
Eq. (12) and the solution of the resulting equa- 
tion for 7, gives 


0.—0; 
1,.= 


tx exp (—@7t rc)2(Rob +r) 

R,+RyetrL—r>d - 
0 86¢*(LR*? ¢?+rR_ +Lr’) 

(14) 

It is necessary to make the use of Eq. (14) less 

laborious. The first step in this direction is to 
write Eq. (10) in dimensionless form as follows 


(rL/R,)=L¢ tan (L¢@). 
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Fic. 2. Values of z for exterior wall 
absorbance equation. 
Total resistance of wall exclu- 
R _ sive of both surface resistances 
R, Surface resistance of 
outdoor side of wall 


In the above equation let the dimensionless 
quantity 


Lo= \ 2 
and also let 


where 


R=the total thermal resistance per unit area of 
either exterior or interior surface of the 
wall, exclusive of both surface resistances. 


The equation can now be written 
(R/ Ry) =z tan Vz. (15) 


In order to make it possible to use Eq. (15), Eq. 
(14) must be written in terms of z and R. In 
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addition, let 
cL=C t=T 
06,—9;=1 =<, 
where 


C=the total thermal capacitance of the wall per 
unit area of either interior or exterior 
surface. 


7 =the time in which the temperature of the air 
in the building must be raised from its 
initial temperature 0; to a comfortable 
temperature 0,. 


a=the heat current that must be supplied per 
unit area of interior wall surface for each 
degree that it is desired to raise the tem- 
perature of the air. 
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When all the above changes have been made, Eq. (14) becomes 


R,+R+R-> 


2RR.+(2R?/z) 


1 


-. (16) 


° (sR? +RR +R?) exp (+27T/RC) 


The determination of z can be made very easy 
by plotting a family of curves having z as ordi- 
nates and R/ R, as abscissas. In order to provide 
ample accuracy, the chart actually used is 18 
inches by 24 inches divided into 1-millimeter 
squares. A simplified reproduction of it is given 
in Fig. 2. For some walls one term of the sum- 
mation is sufficient and five are always ample. 
The quantity represented by the symbol “a” 
should have a name for the sake of convenience 
and, since it is to be used in much the same 
manner as conductance, it will hereafter be 
referred to as the “absorbance” of the wall in 
question. 

Equation (16) is in a sufficiently simple form 
so that the calculation of tables of absorbances 
can be safely left to groups of undergraduate 
students provided those who are reasonably care- 
ful are selected. Once tables of absorbances have 
been prepared, the heat requirements of an 
exterior wall can be found very quickly. For 
instance, if the outdoor temperature is 0°F, the 
indoor temperature 40°F, and it is desired to 


+ 


raise the indoor temperature to 70°F, it is only 
necessary to multiply the conductance of the 
wall by 40, then multiply its absorbance by 30, 
add the two results together and multiply by the 
area of the wall. 

As soon as it can be prepared a paper will be 
published describing a method by which the 
absorbances of more complicated walls such as 
plastered hollow tile with brick veneer can be 
determined with sufficient accuracy to be entirely 
satisfactory for use in connection with the deter- 
mination of required heating plant capacities. 
Another paper will be devoted to the absorbances 
of interior walls and furniture. 

This paper is an abbreviation of a portion of 
the material presented as a doctor's thesis 
at the University of Texas. Especial thanks are 
due to Dr. S. L. Brown for his kindness in 
agreeing to supervise this thesis under very 
difficult circumstances and to Dr. E. G. Keller, 
who is primarily responsible for what knowledge 
the writer has of the kind of mathematics used 
in this paper. 





The Heat Requirements of Simple Intermittently Heated Interior Walls 
and Furniture 


: ELMER G. SMITH 
Agricultural and Mechanical College of Texas, College Station, Texas 
(Received December 17, 1940) 


This paper does for simple interior walls and furniture what the previous paper did for simple 
exterior walls. Unidirectional heat flow is assumed and the boundary conditions are simpler 
than for the simple exterior wall so that the solution is easily expressed in the form of a formula. 
The treatment of the interior surface resistance as a constant quantity is justified. 


N the previous paper, the writer derived an 
equation by which it is possible to compute 
tables that engineers can use to calculate quickly 
the heat requirements of simple intermittently 
heated exterior walls. It is the purpose of the 
present paper to do the same thing for simple 
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interior walls and simple furniture such as church 
pews, theater seats, and table tops. 

Consider a simple interior wall that is initially 
at a uniform temperature throughout, and arbi- 
trarily choose the temperature scale so as to 
make that temperature zero. Also choose the 
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space coordinate so that x=0 at the surface of 
the wall. Thus, as indicated by Fig. 1, the origin 
of temperature is at the initial temperature of 
the interior of the building and the origin of space 
is at the surface of the wall. Also with reference 
to Fig. 1, let: 


L=half the thickness of the wall. 

x =distance from the surface of the wall (x is to be positive 
inward toward the mid-plane of the wall as indicated 
in the figure). 

i=the heat current per unit area at any point in the wall. 

=the temperature, 
in the wall. 

r =the thermal resistivity of the wall. 

c=the thermal capacitance per unit volume of the wall. 


on the arbitrary scale, at any point 


t=time that has elapsed since appreciable heat began to 
enter the room or auditorium. 


the 


As in the case of the simple exterior wall, 
differential equations are 


— (00 /dx)dx = (rdx)1 
— (d1/dx)dx = (cdx)d0/dl 


and the solutions are 


O=A,exp 7 (rcp)'x ]+A2 exp [+(rcp)*x ] 
i=(cp/r)'A, exp [—(rcp)'x ] 
—(cp/r)'A» exp [—(rcp)'x ], 


where p=0/ 01. 
A, and Az:=quantities independent of x. 


As indicated in Fig. 1, the boundary conditions 
are 

when 
when 


i=0 
1=1,1 


x= lL, 
x=), 


where 1=Heaviside’s unit function. 


=the constant heat current that begins 
to enter the wall at time ¢=0 and con- 
tinues thereafter. 


The solution of the equations subject to the 
above boundary conditions gives for the tempera- 
ture at any point in the wall 


t ( L —Z ) . L 
o=ri] — 4] 
rcL 2L 6 


2ri,L, n=+% (—1)" 
_ - > ; exp (—n?2*t/rcL”) 


Tw? n=+1 n? 
x<cos ((L—x/L)nm). (1) 
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Since our immediate objective is to obtain an 
expression for the temperature of the wall surface 
when the air in the building first reaches a 
comfortable temperature we shall let 


T =the time that elapses between the instant 
when appreciable heat first begins to enter 
the room and the instant when the air first 
reaches a comfortable temperature, 


x=, 
0=0,-—0,, 


where 


©, =the temperature of the wall surface on the 
ordinary centigrade or Fahrenheit scale 
when the air in the building first reaches 
the desired temperature. 

(©; =the initial temperature of the building and 
contents also expressed in C®° or F°. 


Then Eq. (1) becomes 


T L 
,—O;= ri + 
te 2 


2ri,L n= te exp (— nial, cL) 
— - } —————— - (2) 


”» ” 


n> n=+ n* 


As in the case of the exterior wall our ultimate 
object is to find an expression for 7, in terms of 7, 
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the constants of the wall, ©; and ©, where 
,=the desired air temperature. 


In order to eliminate ©,, the temperature of the 
wall surface, we may apply Ohm's law to the 
interface between the air and the wall at the 
instant when the air first becomes comfortable 
which gives 


©,,— ©, =1.R,, (3) 
©), 
i.= 
T Ly 2rL 
R.+r + — 
reL 3 a 


In compiling a set of tables, convenience to the 
user should receive careful consideration. In the 
case of the tables for which this paper lays the 
foundation, it will probably happen more or less 
frequently that the exact wall to be heated will 
not be listed in the tables. In order to make it 
easy to decide what wall or walls listed are most 
nearly equivalent to a given wall, it is highly 
desirable to list the total thermal capacitance per 
square foot of surface and the total resistance, or 
conductance, per square foot of surface for each 
wall. Consequently labor will be saved by putting 
Eq. (4) in terms of these quantities. Also, as in 
the case of the previous paper, we must tabulate 
not 1, but 1,/0,—0; which as before will be 
called the absorbance of the wall. Consequently, 
Eq. (4) can be written 


1 
a= -, (5) 
R T  n=« (2R/x*) 
R.+—+—- I aes 
3 C1 nm exp [+(2?/RC)Tn*] 
where 
a=the absorbance of the wall=the number of heat 


units that must be supplied to each square foot of the 
wall for each degree that it is desired to raise the tempera- 
ture of the air in the room or auditorium. 

R=rL=the of the wall between 
mid-plane and its surface per square foot of surface. 


total resistance its 
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where 
R,=the surface resistance of the wall. 


As will be explained later, the treating of the wall 
air interface as a pure resistance does not mean 
that there is any intention to neglect the effect of 
radiation. 

The substitution of the value of ©, from Eq. 
(3) into Eq. (2) and the solution of the resulting 
equation for 7, gives: 


eee ©); 
(4) 
n=+x exp (—n?r*T rcL*) 
n=+1 n- 
C=cL=the total thermal capacitance of the wall 


between its mid-plane and its surface per square foot of 
surface. 


Equation (5) is written in the form that has been 
found most convenient 
calculations. 

It remains to justify the use of the constant R, 
to designate the entire resistance of the wall 
surface to the penetration of heat even including 
radiation from various sources, the temperatures 
of which usually vary at different rates. To 
determine whether R, could be considered a 
constant or would have to treated as a 
function of 7, a considerable number of tests was 
run on a plastered experimental wall 12 inches 
thick, both when bare and when acoustically 
treated. As a result of these experiments and 
some supplementary calculations it was found 
that, for ordinary methods of heating where an 
appreciable part of the heat reaches the walls in 
the form of radiation, R, could be taken as 0.63 
square foot hour degree F per B.t.u., for any 
value of T between } hour and 8 hours. A more 
complete summary 


for making actual 


be 


of these experiments and 
their results will be found in the thesis mentioned 
in the first paper. 

In conclusion, thanks are again due to Drs 
S. L. Brown and E. G. Keller of the University of 
Texas. 
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Effects of Heat Treatment on Field Emission from Metals* 


J. H. DANIEL 
George Eastman Research Laboratories, Massachusetts Institute of Technology, Cambridge, Massachusetts 


(Received April 25, 1941) 


The cylindrical electron projector tube affords a con- 
venient means of investigating the effects of heat treatment 
on field emission from metals. It is found that large, 


unstable field currents, designated as ‘‘anomalous”’ field 


emission, are due to the presence of surface electro-positive 
impurities. As evidenced by the ‘‘blackout”’ of the twinkling 
projector tube patterns, these impurities may be removed 
from tungsten, tantalum, or molybdenum by temperatures 
of 2200°K or greater, but it was found impossible to remove 
them from nickel. While the effect of gas adsorption on 
anomalous field emission is relatively minor, tube pressures 
even as low as 1077 mm Hg increase the frequency of 
“‘breakdown.”’ This phenomenon is probably caused by 


INTRODUCTION 


IELD emission of electrons from metals 

occurs when the application of a strong 
electrostatic field so lowers and thins the 
potential barrier at the surface of the metal 
that electrons escape through the barrier by 
virtue of wave-mechanical penetration. One of 
the earliest recognized characteristics of field 
emission currents was their extreme dependence 
on the heat conditioning of the emitting surface, 
particularly during the initial raising to high 
temperatures. 

Millikan and Eyring! made extensive measure- 
ments of current as a function of voltage for 
various amounts of heat treatment on thoriated 
tungsten filaments.2 They found that a semi- 
logarithmic plot of current vs. voltage for a raw 
wire with no previous heating showed discon- 
tinuities in a decreasing slope as the voltage was 
raised to a maximum value. As the voltage was 
reduced, a comparatively smooth curve with 
greater average slope for the same voltages (or 
the currents) was followed. This 
smoother curve was in general reversibly repro- 


even same 


* Part of a thesis presented for the degree of Doctor of 


Science, from the Department of Physics, Massachusetts 
Institute of Technology, October 1940. 

1 R.A. Millikan and H. Eyring, Phys. Rev. 27, 51 (1926). 

2 No mention is made of any effect which might be 
expected due to activation of the filaments. If none was 
observed, the fact would stand in agreement with the 
experiments of A. J. Ahearn, Phys. Rev. 50, 238 (1936). 
A probable explanation is given in reference 15. 
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rearrangement or disruption of impurity layers, or even of 
the metal itself, due to positive ion bombardment. Making 
surfaces microscopically smooth greatly decreases the 
magnitude of anomalous field emission currents, and 
therefore the likelihood of breakdown, but cannot take 
the place of sufficiently severe heat treatment in removing 
the anomalous type of emission completely. In the course 
of these experiments, calculations were made which make 
possible the prediction of the tension on the projection 
tube filament necessary to prevent violent vibration of the 
filament as high voltages are applied to the tube. Since 
such cylindrical tube geometry is of common use in experi- 
mental work, the results of these calculations are given. 


ducible if the first maximum values of current 
and voltage were not exceeded. After a period 
of heating, this same type of behavior was again 
exhibited, but for temperatures above 1100°K, 
progressively higher voltages were required to 
produce the same currents. 

From tables presented by Millikan and 
Eyring,! it is evident that the voltage necessary 
to produce a given current, defined as the 
“critical voltage,” increases very markedly for 
heat treatments at or above about 1900°K. 
They explain this as being due to removal of 
surface impurities, plus a rounding off of the 
surface peaks from which most of the emission 
took place (due to concentration of the electro- 
static field around these peaks, and evidenced by 
fluorescent spots on the glass walls of the tube) 
both by the heating and by the more energetic 
positive ion bombardment at higher voltages. 
The drastic increase in critical voltage was 
apparently not permanent, as the critical voltage 
decreased again with long standing. Several 
cases in which the current suddenly increased 
discontinuously by factors of 100 or 1000 were 
observed. 

Observations of this kind, while having no 
significance to fundamental theory of metals, 
are of practical importance in determining the 
kind of metal and ‘“‘conditioning”’ process which 
should be given high voltage tubes to minimize 
the possibility of breakdown due to field emis- 
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sion. Measurements such as Millikan made are 
extremely time consuming and ‘tedious,” but the 
electron projector tube offers a quicker way for 
obtaining the desired information. 


EXPERIMENTAL APPARATUS 


The type of projector tube used is shown in 
Fig. 1. The calcium 
tungstate and willemite sprayed on the glass 
wall and covered with a potassium silicate 
coating (to prevent bombardment of the filament 
with particles detached from the screen) as 
Martin.’ The potential 
closely follows the collector potential up to 


fluorescent screen is of 


described by screen 
perhaps 10,000 volts by virtue of its secondary 
emission properties. Emission of electrons from 
the wire takes place radially due to the geometry 
of the electrostatic field and its concentration 
around the filament. Thus the electrons striking 
the screen give an emission image of the wire sur- 
face which is greatly magnified in circumference.® 

For the high voltages necessary to obtain field 
emission patterns, it was found that quite often 
the filament suddenly broke into violent vibra- 
tion, making observations impossible. Since this 
tube geometry is of common use in experimental 
work, and the same difficulty will limit the 
range of any measurements made with voltages 
upward of 1000 volts, some time was spent in 
investigating the phenomenon, and information 
was obtained which enabled the behavior of a 
given tube to be predicted. 

Experimental facts easily observed were: As 
the voltage on a particular tube was raised, a 
reproducible, critical value was reached at which 
the filament quite abruptly broke into fairly 
large amplitudes of vibration. The frequency 

4S. T. Martin, Phys. Rev. 56, 947 (1939). 

* As given by the equation 

F= V 
(r+x) In R/r 


where F is the potential gradient or field at distance x 
measured radially outward from the wire surface, V the 
potential difference between wire and screen, and r and R 
the respective radii of the wire and screen. 

5 Distortions at the wire surface (where the field is 
greatest) of the cylindrical symmetry of the field, due to 
its local concentration about the small projections from 
which field emission takes place, account in large measure 
for the fact that field emission images consist of elongated 
bright spots rather than very thin lines around the screen 
circumference. See Fig. 2. 
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Fic. 1. Cylindrical electron projector tube. 


was not necessarily 60 cycles per second nor any 
multiple thereof; the vibration did not take 
place in a fixed plane, but in one which rotated 
in azimuth around the stationary filament axis; 
almost never was any but the fundamental mode 
present. Such characteristics seemed to point to 
a condition of unstable equilibrium being sud- 
denly attained by the stationary wire as the 
force on it produced by the electrostatic field 
exceeded a critical value. With this in mind, 
the forces on the wire were computed,® and the 
following expression obtained for the voltage 
necessary to initiate vibration: 


R 
V>6.8X104(7/1)R logy (-). (1) 


r 


where V is the voltage between screen and 
filament, 7 the tension in grams on the filament, 
l the length of the filament in centimeters, R 
the radius of the screen in centimeters, and r 


® J. H. Daniel, ‘Studies of field emission,’’ Sc. D. Thesis, 
Massachusetts Institute of Technology (1940), Appendix I. 
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the radius of the filament in centimeters. 
Solving for the tension necessary to prevent 
vibration: 
VP 
T> 2.16 X 10-"—_______ (2) 
R*(logio (R r) ?? 


For a wire not accurately centered, the effect® 
on (1) and (2) is small and tends to lessen the 
tension necessary to prevent vibration. 

Using (1) or (2) in conjunction with Blodgett 
and Langmuir’s’ specifications for design of a 
helical tungsten spring® to give maximum safe 
tension at temperatures up to 2800°K, where 
the tension vanishes (due to expansion of the 
filament), the onset of vibration can be predicted 
for any temperature. If the filament is not to be 
heated to 2800°K, greater tension may be applied 
to hold off vibration to higher voltages. If high 
voltage is to be applied only when the filament 
is cold or at low temperatures, a short stiff 
spring may be used which goes slack at a safe 
temperature, but exerts a strong tension at lower 
temperatures. The behavior of all tubes con- 
structed carefully with hydrogen-fired springs’ 
according to the above specifications was found 
to confirm the theory from which (1) and (2) 
were derived. 

To extend the specifications to filament metals 
other than tungsten, it was found sufficiently 
good to calculate the allowable stretch of the 
tungsten spring as the fraction of that required 
for a tungsten filament represented by the 
product of the ratio of the melting point of the 
metal to that of tungsten by the ratio of their 
respective linear coefficients of expansion. This 
gives the following approximate stretches: tanta- 
lum stretch = 0.83 tungsten stretch, molybdenum 
stretch=0.77 tungsten stretch, nickel stretch 
= 1.06 tungsten stretch. 

All filaments used were of 0.001 inch diameter 
except those of molybdenum, which were of 
().002 inch diameter. Temperatures of the tung- 
sten wires were computed from the diameter 
(determined by accurate weighing of a measured 
length) and the current input, using the Jones- 


7K. B. Blodgett and I. Langmuir, Rev. Sci. Inst. 5, 321 
(1934). 

§ Specifications for molybdenum springs may be found 
in an article by G. L. Tawney, Rev. Sci. Inst. 10, 152 
(1939). 
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Langmuir tables.* Temperatures of the tantalum 
wires were computed similarly using tabulated 
values of A’ in Table I of Malter and Langmuir.'” 
Temperatures of wires of other metals were 
measured with a special micropyrometer con- 
structed for use with filaments of small diam- 
eter," and calibrated on a lamp standardized by 
the National Bureau of Standards. Measure- 
ments of temperature, made on the tungsten 
and tantalum wires with this pyrometer and 
compared with the calculated values, indicated 
a probable discrepancy of less than +20°K for 
values taken from a graph of individual measure- 
ments, and less than +50°K for the individual 
measurements themselves. Filaments of about 
20 centimeters length were used, so corrections 
for end losses due to heat conduction through 
the leads were negligible. For the type of 
observations now to be described, such tempera- 
ture measurements were than 
enough. 


more accurate 


OBSERVATIONS 


Figure 2 shows typical field emission patterns 
obtained from a raw tungsten General Electric 
218 wire (No. 1 in Table I) at various stages of 
heat conditioning. In each case the filament 
was kept at a fixed temperature, usually for 
ten minutes, then brought to room temperature. 
The voltage necessary to produce the pattern 
was applied immediately, and the pattern 
photographed. For temperatures up to about 
1600°K the pattern was far from being the 
static representation necessarily shown in a 
photograph, but consisted of numerous brilliant, 
twinkling spots which continually appeared, 
disappeared, and reappeared. This activity was 
quite intense for a period of several minutes, 
but if the field remained on, it gradually di- 
minished. The pattern meanwhile became less 
brilliant, but- never became completely stable 
even after several hours. If the voltage was 
raised and then lowered, patterns obtained 
during the decreasing range of voltages were 
less brilliant and probably somewhat less active 


*H. A. Jones and I. Langmuir, Gen. Elec. Rev. 30, 310, 
354, 408 (1927). 

1 L. Malter and D. B. Langmuir, Phys. Rev. 55, 743 
(1939). 

1M. L. Benson, Master’s Thesis, Massachusetts Insti- 
tute of Technology (1935). 
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Fic. 2. Electron projection tube patterns produced by 
anomalous field emission from a 0.001-inch diameter raw 
polycrystalline General Electric 218 tungsten wire at room 
temperature. 

I. After heat treatment incidental to baking of tube at 800°K 
(pattern is similar if tube is not baked at all). Fluorescent screen 
potential =5000 volts. 

Il. After 10 minutes at 1200°K. Screen potential =5000 volts. 

Ill. After 10 minutes at 1500°K. Screen potential =5000 volts. 

IV. After 10 minutes at 1700°K. Screen potential = 7000 volts. 

V. After 10 minutes at 1800°K. Screen potential =9000 volts. 

VI. After 10 minutes at 1900°K. Screen potential = 10,000 volts. 


(by activity is meant the amount of flickering 
present) than those for the same voltage values 
during the initial voltage increase. This corre- 
sponds qualitatively to the current-voltage 
curves obtained by Millikan and Eyring, and 
the process was called by them “current 
conditioning.”’ 

For temperatures greater than 1600°K, the 
pattern became progressively less brilliant (or 
-higher voltages were required to produce a 
pattern), and the emission spots more stable and 
less numerous. Not infrequently as the voltage 
was raised, an extremely brilliant spot would 
suddenly appear, gradually diminish in brilliance 
for a minute or two, but remain permanently in 
the pattern. Such spots must have indicated a 
disruption of the surface, or of the impurity 
layers on the surface”; they were sometimes so 


2 Such disruptions of the tungsten surface were actually 
observed with the optical microscope in the case of field 
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violent as to cause a sort of ‘‘recoil’’ of the 
filament. They correspond to the breakdown 
phenomena mentioned in the literature. 

The most important characteristic of this heat 
conditioning, however, is the lack of pattern, or 
complete “blackout” obtained after heating to 
2000°K, for voltages several times as large as 
those previously required to produce a pattern. 
For the wire of Fig. 2, 20,000 volts were 
insufficient to produce a single trace of light from 
the screen. If higher voltages could be applied 
without vibration setting in, a breakdown usually 
occurred, which, in most cases, healed up after 
another ten minutes heating at 2000°K. Heating 
periods of one minute or less at 2000°K were, 
in general, insufficient to produce blackout, but 
if the temperature was elevated to 2300°K or 
2400°KK, such short periods were sufficient. 
Periods of several hours below 1900°K did not 
produce blackout. The emission giving rise to 
the twinkling patterns existing before the occur- 
rence of blackout will be called ‘‘anomolous’”’ 
field emission, to distinguish it from the steadier, 
more predictable emission which took place 
after blackout, but was so small as to be insuffi- 
cient to cause any visible fluorescence. 

Figure 3 shows the heat conditioning history 
of a 0.002” Callite Type 400 RI molybdenum 
wire (No. 9). At 1810°K a blackout was obtained, 
but on raising the temperature further to 1850°K, 
spots again appeared. These remained until the 
2260°K heating period (enough evaporation took 
place at this temperature to darken the screen 
slightly), whereupon blackout appeared com- 
plete. This early blackout, followed by subse- 
quent reappearance of anomalous field emission, 
occurred at 1750°K in one (No. 10) of the other 
two molybdenum wires (No. 5 and No. 10) 
tested. The other (No. 5) exhibited permanent 
blackout at 1900°K. The accompanying thermi- 
onic pictures (with subscripts A) were taken 
immediately following the field emission pattern 
pictures. There was usually a fair correspondence 
in position of bright bands in thermionic patterns 
emission from sharp points. Positive ion bombardment, 
concentrated on the small emitting areas by the electro- 
static field geometry, could cause such intense local heating 
as to reduce the tensile strength of the metal below the 
strong electrostatic forces to which it is subjected. 

‘8 That the corresponding surface field existed even 


though the fluorescent screen perhaps did not rise above 
10,000 volts was assured by the nickel anode. 
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to positions of spots in field emission. patterns. 
These bands also disappeared when final black- 
out occurred for anomalous field emission (they 
did not disappear for the ‘‘premature”’ blackout 
which was observed only for molybdenum), and 
a uniform pattern with longitudinal die marks 
(which did not appear with wires polished 
microscopically smooth) as the chief feature 
remained. 

After blackout had been produced, admitting 
air to the tube containing wire No. 9 and then 
pumping it down again to a pressure of 10~® or 
10-7 mm Hg without even baking it failed to 
remove the blackout. Even after the tube lay 
around in the laboratory at atmospheric pressure 
with a cork in the tubulation for a week, no 
anomalous field emission appeared when it was 
again evacuated. When it had remained at 
atmospheric pressure on the vacuum system 
for another week, the tube showed small spots 
as in Fig. 3-VI upon evacuation and application 
of 10,000 volts. An idea of the dimness of these 
spots may be obtained from the fact that this 
was a 60-second exposure as compared to the 
usual one- to ten-second exposures. Heating the 
filament to only 1200°K removed the spots. 








Fic. 4. Electron projector tube patterns produced by 
anomalous field emission from a 0.001-inch diameter raw 
polycrystalline Driver-Harris Grade ‘‘A”’ nickel wire at 
room temperature. 


I. No heat treatment. Screen potential =3000 volts. 
II. After 10 minutes at 1530°K. Screen potential =5000 volts. 


After two more weeks at atmospheric pressure 
on the vacuum system, one spot appeared, but 
was removed by a temperature of 1100°K. 
Molybdenum wire No. 10 exhibited the same 
sort of behavior, as did tungsten wire No. 2, 
slight thermionic bands accompanying the field 
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Fic. 3. Electron projector tube 
patterns produced by anomalous 
field emission (by thermionic emis- 
sion in pictures with subscripts A) 
from a 0.002 inch diameter raw 
polycrystalline Callite type 400 RI 
molybdenum wire at room temper- 
ature. 


I. After 10 minutes at 1200°K. Screen 
potential =4000 volts. 

II. After 10 minutes at 1700°K. Screen 
potential =6000 volts. 

I-A. Thermionic pattern for same con- 
ditions as in Fig. 3-II. Screen potential = 
4000 volts. Temperature =1700°K. 

III-A. Thermionic pattern after 10 min- 
utes at 1810°K. Screen potential =4000 
volts. Temperature = 1800°K, 

IV. After 10 minutes at 1900°K. Screen 
potential = 11,000 volts. 

IV-A. Thermionic pattern for Fig. 3-1V. 
Screen potential = 4000 volts. Temperature 
=1900°K. 

V-A. Thermionic pattern after 10 minutes 
at 2200°K. Screen potential =4000 volts. 
Temperature =1900°K. Note longitudinal 
die marks. 

VI. After two weeks at atmospheric pres- 
sure. Screen potential = 10,000 volts. 





VOLUME 12,"AUGUST, 1941 


649 














I sn 


Il 





Fic. 5. Electron projector tube patterns produced by 
anomalous field emission from a 0.001-inch diameter raw 
polycrystalline Driver-Harris Grade ‘‘A”’ nickel wire at 
room temperature. 


No heat treatment. Screen potential =4000 volts. 


I. After 10 minutes at 1200°K. Screen potential =4000 volts. 
Il. After two hours at 1600°K. Screen potentia! =5000 volts. 


emission spots. Wire No. 2 was then removed 
from the tube and placed in glass tubing. After 
four been 
handled only with cotton-gloved hands. Upon 


days it was remounted, having 
evacuation it behaved exactly as though under- 
going heat treatment for the first time. Two 
other tungsten wires, one tantalum, and one 
molybdenum wire behaved in the same manner. 

Figures 4-I and 4-II show the patterns of a 
0.001". Driver-Harris Grade ‘A’”’ nickel wire 
after temperature and at 
1530°K, respectively. The wire (No. 3) broke 
while cold under 5000 volts. Figure 5 shows the 


treatment at room 


history of a second nickel wire which broke at 
1650°K. 

Figure 6 shows the patterns (badly under- 
exposed) obtained from a tungsten wire polished 
microscopically smooth by a mechanical grinding 
described by Nelson.“ The patterns 
were very dim, and consisted of only a few spots, 
but the blackout was quite definite, and occurred 
at approximately the same temperature (1900°K) 
as for the unpolished tungsten wires. On the 


process 


smooth wire the spots were much more stable, 
and the 
noticeable. 
Figure 7 shows tantalum patterns. With the 
two tantalum wires (Nos. 8, 11) investigated a 
blackout also occurred, although at a slightly 
higher temperature—about 2300°K. However, 


effects of current conditioning less 


4R. B. Nelson, 
filaments,” Ph.D. 
Technology (1938). 


structure of 
Massachusetts 


“The grain 


tungsten 
I hesis, 


Institute of 
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the blackout was not quite so complete as with 
tungsten and molybdenum, for spots could be 
obtained from both wires by applying sufficiently 
high voltages (from 15,000 to 20,000 volts) even 
after ten minutes at 2600°K. Both wires broke 
at 2800°K. 

Table I further the results 
obtained from tungsten, tantalum, molybdenum, 
and nickel wires. 


summarizes 


CONCLUSIONS 


The less spectacular anomalous field emission 
pattern given by the microscopically smooth 
tungsten wire, and its disappearance at 1900°K, 
show that although surface roughness has a 
great deal to do with the magnitude of anomalous 
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Fic. 6. Electron projector tube patterns produced by 
anomalous field emission from 0.001-inch diameter G. E. 
218 tungsten wire polished microscopically smooth. 











. No heat treatment. Screen potential =5000 volts. 

I. After 10 minutes at 1650°K. Screen potential =6000 volts. 

II. After 10 minutes at 1850°K. Screen potential =9000 volts. 
II-A. Thermionic pattern for Fig. 6-III. Screen potential =4000 
volts. Temperature =1800°K. 


I 
I 
I 
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TABLE I. 


Wiki AvG. PRESSURE BLACKOUT MELTING 
No. DESCRIPTION DuRING TES! TEMPERATURE POINT REMARKS 
1 0.001-inch 4x1075 2000°K for 3640°K 
G. E. 218 10 minutes 
tungsten 
2 Same as No. 1 310-8 mm 1870°K for 3640°K 
Hg 10 minutes 
4 Same as No. 1, 21075 1900°K for 3640°K 
but polished 10 minutes 


microscopically 
smooth 


8 0.0011-inch 2%16"* 2400°K for 3120°K Blackout not so definite nor so 
Fansteel Ta 10 minutes complete. Wire broke under ten- 
sion at 2800°K. 
11 Same as No. 8 not measured 2200°K for 3120°K Same as for No. 8. 
10 minutes 
5 0.002-inch 2x 10-8 1900°K for 2890°K 
Callite Type 10 minutes 
400 RI Mo 
9 Same as No. 5 not measured 2260°K for 2890°K Premature blackout occurred at 
10 minutes 1810°K. 
10 Same as No. 5 not measured 2200°K for 2890°KK Premature blackout occurred at 
10 minutes 1750°K. 
3 0.001-inch 6x 1077 1725°K Anomalous field emission - still 
Driver-Harris definitely present after 10 min. at 
Grade ‘*A”’ 1530°K. Wire broke while cold 
nickel under 5000 volts. 
6 Same as No. 3 2107-5 1725°K Anomalous field emission still de- 


finitely present after 2 hours at 
1600°K. Wire broke while sagging 
at 1650°K. 
7 Same as No. 3 3X 1075 1725°K Wire broke while cold under 5000 
volts (defective weld?). 


~ 


Note.—Wires were heated at each new temperature (usually 100°K higher) for at least ten minutes—longer if necessary to bring pressure back 
to former value. Pressures were measured with a tantalum plate ion gauge calibrated on air against a McLeod gauge in the range 10-* to 10-* mm 
Hg, and checked for Barkhausen oscillations and photoelectric emission from the plate. All observations and measurements were made with the 
tubes on the vacuum system and liquid nitrogen around the trap. 


field emission currents, it is not a factor in anomalous field currents, but is no more im- 
determining the blackout temperature. Thus a_ portant than sufficiently severe heating of the 
smooth surface is advantageous in lessening surface, which alone can remove such currents 
completely. 

The persistence of blackout after air has been 
admitted to the emitting surface shows that 
anomalous field emission is not due to adsorbed 
gases. This is a fact which could not be assumed 
without such proof since the blackout tempera- 
tures are in the neighborhood of those required 
to remove some of these gases from metallic 
surfaces. The logical cause is thus narrowed to 
electro-positive surface impurities‘such as the 
alkalis, which tend to lower the work function, 
and which normally require some time to reach 
the filament unless it is handled or brought into 
contact with objects such as other metals or 
glass (e.g., tungsten wire No. 2). Further 
investigations of surface conditions using point 








Fic. 7. Electron projector tube patterns produced by 
anomalous field emission from 0.001-inch diameter Fansteel : - . : 
tantalum wire at room temperature. projector tubes'® make it apparent that while 


I. No heat treatment. Screen potential =3000 volts. . ‘ . x 
II. After 10 minutes at 2000°K. Screen potential =10,000 volts. 4% J. H. Daniel, Phys. Rev. (to be published). 
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surface electro-positive material is necessary 


and sufficient for the production of anomalous 
field emission, with any but the best surface 
conditions'® such material is accompanied by 
the presence of adsorbed gas. 

In contrast to the lack of dependence of 
blackout (or removal of electro-positive ma- 
terials) on vacuum conditions, the phenomenon 
of breakdown certainly occurs more frequently 
under relatively poor vacuum conditions (pres- 
sures of 10-7 mm Hg or greater). This can be 
explained if local disruption of the surface due 
to violent positive ion bombardment is assumed™ 
to be the primary cause of breakdown. The 
same assumption explains the greater frequency 
of breakdown observed with surfaces not sub- 
jected to heat treatment, since the greater 
emission from the impurity layers produces more 
positive ions. That an actual disruption of the 


It is also demonstrated strikingly that total pressure 
is not always a significant criterion of surface condition. 


surface takes place, and in some _ instances 
extreme local heating occurs, has been shown 
with field emission from points.'® With such 
points it was proved also that the production 
of blackout by heating is not due to smoothing 
of the surface by surface tension (since this 
occurs at much higher temperatures)—a fact 
which is indicated here by the return of anoma- 
lous field emission patterns upon contamination 
after blackout has been produced. 

The metals tungsten, tantalum and molybde- 
num can be rid of anomalous field emission by 
short heating periods at 2200°K or above. It is 
impossible to clean nickel sufficiently at any 
temperature below its melting point to get rid 
of anomalous field emission. 

The author wishes to express his appreciation 
to Professor W. B. Nottingham for suggesting 
this problem in field emission. Thanks are also 
due Mr. L. W. Ryan for aid in glass blowing 
operations. 





Here and There 








Summer Course in Photoelasticity 


To meet the growing interest in photoelasticity, the 
Department of Mechanical Engineering of the Massa- 
chusetts Institute of Technology, Cambridge, Massa- 
chusetts, is sponsoring a short summer course dealing with 
this method of stress analysis. This course will extend from 
July 28 to August 8. Its inclusion in the summer session is 
primarily for the benefit of those engaged in industry or 
other occupations which prevent attendance at the In- 
stitute’s regular winter session. 

The thirteenth semi-annual meeting of the Eastern 
Photoelasticity Conference was held at the Massachusetts 
Institute of Technology, under the auspices of the Depart- 
ment of Mechanical Engineering, on June 12, 13, and 14. 


* 


Applied Mechanics Courses Offered 


A staff of sixteen experts from Europe and America, 
representing ‘“‘the best men available in the fields of 
mathematics, engineering and physics,”’ has been selected 
to take charge of summer instruction and research in 
applied mechanics at Brown University, - Providence, 
Rhode Island, according to an announcement made by the 
University’s news bureau. The summer session, which 
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began on June 23 and is to continue through September 13, 
is designed to help overcome a bottleneck in the supply of 
men with highly specialized knowledge who are needed in 
defense industry. Dean Roland G. D. Richardson of 
Brown’s Graduate School is the Director of the program. 


* 


Chicago University’s Fiftieth Anniversary 


More than 160 of the nation’s leading scientists and 
scholars, including thirty-two distinguished men and 
women who will be awarded honorary degrees, will report 
basic achievements and advances in learning in a five-day 
series of symposia sponsored by the University of Chicago. 
The symposia will be held on the Midway Quadrangles in 
the week climaxing the celebration of the University’s 
Fiftieth Anniversary, beginning September 22, according 
to a recent announcement. Thirty-nine universities, in- 
cluding six in foreign nations, and fifteen museums, research 
organizations and government agencies will be represented 
in the symposia. They will deal with newest fundamental 
advances in the biological, physical, and social sciences, 
the humanities, law, business, religion, and social service, 
in keeping with the theme of the University’s celebration, 
“New Frontiers in Education and Research.’’ Sessions in 
the biological and physical sciences will be held jointly 
with the American Association for the Advancement of 
Science. Dr. Frederic Woodward, vice president emeritus, 
is Director of the Fiftieth Anniversary. Robert A. Millikan 
and Ernest O. Lawrence, Nobel prize winners in physics, 
are among the speakers. 
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New Society Holds First Meeting 


The American Society for X-Ray and Electron Dif- 
fraction starts its existence with a charter membership of 
124. The officers elected for 1941 are: President, M. L. 
Huggins of the Eastman Kodak Company; Vice President, 
B. E. Warren of the Massachusetts Institute of Technology; 
Secretary-Treasurer, George Tunell of the Geophysical 
Laboratory, Washington, D. C. 

The first meeting of the Society was held at Gibson 
Island, Maryland, in cooperation with Section C of the 
American Association for the Advancement of Science. The 
program scheduled was as follows: 

Monday, July 28 
10:00 a.m. B.E. WARREN. Theory and Practice of Particle Size 
Determination. 

7:30 p.M. M. J. BUERGER. New Experimental Techniques in Crys- 

tal Structure Analysis. 


Tuesday, July 29 
10:00 a.m. A. L. Patterson. Fourier Series as Applied to Crystal 
Structure Determination. 
7:30 p.m. C.S. Futter. X-Ray Studies of Linear Polymers of 
Known Chemical Constitution. 
Wednesday, July 30 
9:00 A.M. N.R.C. COMMITTEE ON X-RAY AND ELECTRON Dir- 
FRACTION 
10:00 a.m. P. Desye AND P. P. Desye. Discussion of Theoretical 
and Experimental Electron Diffraction Methods. 
2:00 p.m. Business meeting. 
7:30 p.m. E.W. HuGues. Recent X-Ray and Electron Diffraction 
Work at the California Institute of Technology. 
Thursday, July 31 
10:00 a.m. L.O. Brockway. Recent Electron Diffraction Work on 
Gas Molecules. 
2:00 p.m. M.L. Futter. Applications of Electron Diffraction to 
Industrial Problems. 
Friday, August 1 
10:00 a.m. E. F. Burton. Applications of the Electron Microscope. 
2:00 p.m. V. K. Zworykin, J. HILLIER, AND A. W. VANCE. Pre- 
liminary Report on 300,000-Volt Electron Microscope. 


- 
Sloan Fellowship Award 


W. Endres Bahls, who is in charge of development and 
design work in connection with special radio tubes at the 
Harrison, New Jersey, laboratories of the RCA Manu- 
facturing Company, has been awarded an Alfred P. Sloan 
Foundation Fellowship for a year of advanced study of 
industrial problems at the Massachusetts Institute of 
Technology. One of eleven engineers chosen from. all 
United States industry to receive the honor, Mr. Bahls 
will spend one year at the Institute attending the Business 
and Engineering Administration School. 

The purpose of the Fellowship is to permit men actively 
engaged in industry to study such subjects as finance, ac- 
counting, labor relations, economics, law, and marketing, 
stressing the sociological aspects of modern industry. 


* 
Modern Plastics Competition Announced 


Modern Plastics Magazine is now sponsoring its Sixth 
Annual Modern Plastics Competition. The sponsors of this 
annual competition invite all plastic-using firms, designers, 
molders, laminators, fabricators, materials suppliers, 


VOLUME 12, AUGUST, 1941 


machinery and die makers to participate. Any plastic 
object or product will be deemed eligible if it has been 
designed or has reached the market since September 1, 
1940. The competition and concurrent display of entries 
will be held in the main headquarters of Modern Plastics, 
Chanin Building, 122 East 42nd Street, New York. The 
final deadline for all entries is September 8. 


* 


Scientific Editor Receives Appointment 


Samuel Glasstone has been appointed Scientific Editor 
of Princeton University Press, it was announced recently 
by Joseph A. Brandt, Director of the Press. Dr. Glasstone 
was born in London, England, and was educated in the 
University of London, where he obtained the degrees of 
Doctor of Philosophy and Doctor of Science. He taught 
physical chemistry in the University of London, the Uni- 
versity College of the Southwest of England and the 
University of Sheffield. Since 1939 he has been research 
associate in chemistry at Princeton University. Dr. 
Glasstone is internationally known for his researches in 
physical chemistry and particularly in electrochemistry. 

The object of Princeton University Press in appointing 
Dr. Glasstone as its Scientific Editor is to encourage and 
facilitate the publication of scientific books, particularly 
works of a specialized character which may lie outside the 
purview of the ordinary commercial publishers. 


* 


Films on Flying Available 


The first basic films on flying available to educational 
institutions, prepared by Erpi Classroom Films, of Long 
Island, New York, in collaboration with the University of 
Chicago, were released in June. The films were produced 
under the direction of Dr. Reginald J. Stephenson, As- 
sistant Professor of Physics at the University of Chicago, 
and Lieutenant Walter T. Brownell, Director of the 
University’s student pilot training program. The pre- 
sentation is in two reels, ‘Theory of Flight” and ‘Problems 
of Flight.” 

* 


Calendar of Meetings 


August 
27-29 American Institute of Electrical Engineers, Yellowstone Na- 
tional Park 


September 
8-12 American Chemical Society, Atlantic City, New Jersey 
8-10 American Astronomical Society, Williams Bay, Wisconsin 
22-26 American Association for the Advancement of Science, Chicago, 
Illinois 
5-26 Society of Automotive Engineers, Milwaukee, Wisconsin 


N 


October 

6-10 National Safety Council, Chicago, Illinois 

12-15 American Society of Mechanical Engineers, Louisville, Ken- 

tucky 

13-16 Society of Motion Picture Engineers, New York, New York 
15-17. American Society of Civil Engineers, Chicago, Illinois 
20-24 American Society for Metals, Philadelphia, Pennsylvania 
20-24 American Welding Society, Philadelphia, Pennsylvania 
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Letter to the Editor 








Deterioration of Photographic 
Films before Development* 


Wayne T. SpROULL 
General Motors Research Laboratories, Detroit, Michigan 
May 8, 1941 


HOTOGRAPHIC films and plates are commonly sold 
with an expiration date stamped on the package which 
implies a warning that deterioration with time is to be 
expected before development. Discussion of the nature and 
causes of this deterioration is meager in the publications 
on photography. High temperature and high humidity 
accelerate the deterioration, but there may be other factors 
involved because film and plates deteriorate slowly even 
when stored in a cool dry place. 

If cosmic rays were one of the principal factors re- 
sponsible for the deterioration of films and plates, the 
that 
discover a way to prevent such deterioration would be 


prospects photographic manufacturers could ever 
very poor. Hence the negative character of the results of 
these experiments is in a sense encouraging because it 
confirms the generally accepted theory that film and plate 
deterioration before development is due to humidity and 
temperature effects which the manufacturers may some 
day overcome. 

Millikan 


experimental data that the total energy per square centi- 


and Cameron! have computed from their 
meter per second that flows into the top of the earth’s 
atmosphere in the form of cosmic rays is about one-tenth 
of the total energy coming into the earth’s atmosphere in 
the form of starlight and heat. One would expect film or 
plates exposed to starlight for a few months to show signs 
of fogging. Old film when developed behaves in some ways 
like slightly fogged film. It seems reasonable, then, that 
cosmic rays might be partly responsible for the pre-de- 
velopment deterioration of photographic films and plates. 

Upon this theory, one might expect x-ray films to de- 
teriorate somewhat more rapidly than films intended for 
exposure to ordinary light, because their thicker double 
emulsions, designed to make them more sensitive to x-ray 
exposure might also make them more sensitive to cosmic- 
ray exposure. In order to test this theory, some Eastman 
“Verichrome”’ roll film and some Agfa ‘‘Non-screen”’ x-ray 
film was purchased in Detroit on February 3, 1939. One- 
half of each type of film was sent to Mr. Harry S. Donald 
of the Calumet and Hecla Consolidated Copper Co. and 
Allied Companies, of Calumet, Michigan. Mr. Donald 
kindly stored these films on February 27, 1939, 2057 feet 
vertically beneath the surface of the ground in a shaft 
just south of the one used by Wilson? in his cosmic-ray 
experiments. Both shafts are in the same type of rock, 
and Wilson’s data indicate that at this depth, the cosmic- 
ray intensity is only about one twenty-thousandth of that 
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at’ ground level. Meanwhile, the remaining films were 
stored at an ordinary place in Detroit. 

On March 28, 1941, the films were removed from the 
mine and on April 2, 1941, some of these mine-stored 
x-ray films were developed along with some of the x-ray 
films stored in Detroit, and a piece of new Agfa “‘Non- 
screen” film. These films were developed and fixed together. 
Figure 1 shows a print made from the three films. Both the 
film stored in the mine and that stored in Detroit appear 
fogged when compared to the new film, and the fogging 
is about equal in both old films. The conclusion is obvious 
that cosmic-ray exposure is not a major factor in the 
deterioration of x-ray films before development. 

An additional test of the x-ray films was made by radio- 
graphing some steel strips laid in a criss-cross pattern over 
all three types of film simultaneously. This radiograph 
shows that, in addition to the fogging effect of age (causing 
the very lightly exposed portions of the negative to appear 
darker), there is a loss of sensitivity with age, shown by the 
fact that the heavily exposed portions of the old film are 
not as dark as the equally heavily exposed portions of the 
new film. One might suspect that cosmic rays could be 
partly responsible for fogging old film, while the loss of 
sensitivity might be due to other causes. However, the 
radiograph shows that the mine-stored film behaves in all 
respects like old film, which indicates that cosmic rays 
are not a major factor in causing either phase of the 
deterioration. 

On April 16, 1941, tests were made on the Verichrome 
roll film. The film was loaded in a camera which was used 
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Fic. 1. This print was made from the three different films by a 
single exposure. The fact that the mine storage film in the center has 
fogged with age to practically the same degree as the ordinary storage 
film is evident. 


to photograph a small illuminated sign in an otherwise 
dark room. The sign served to mark the film permanently 
so it could be identified as new film or as film stored in the 
mine or in Detroit. The results were in agreement with 
those obtained in the x-ray film experiments. It is evident 
that the mine storage film is inferior to the new film and 
that it behaves like the film stored in Detroit. 

* This work was carried out on the author's own initiative and was 
not connected with the work of the General Motors Laboratory. 


1R. A. Millikan and G. H. Cameron, Nature 121, 26 (1928). 
2 Volney C. Wilson, Phys. Rev. 53, 337 (1938). 
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